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ABSTRACT 


The results of a/ program to study the analysis and design of composite 
materials and structures is reported. Emphasis was placed upon three 
major areas: The definition of design criteria for laminates including 
studies of basic .failure mechanisms; the definition of unique design* concepts 
to enhance the beneficial characteristics. of composite materials and to utilize 
them 'in structures; and the analysis of composite materials property test 
te chniques. 
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FOREWORD 


This document .is the annual report on the program entitled "Study of 
the Relationship of Properties of Composite Materials to Properties of 
The'ir Constituents". The program was performed for the National Aero- 
nautics and Space Administration, under Contract NASw-1377 and was 
monitored by Dr. R. W. Leonard of the NASA Structures Research Divi- 
sion. 
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•INTRODUCTION 


Development of -.methods of analysis of the strength and stiffness of fib- 
rous composites has advanced to the point where it is feasible to establish 
rational initial design procedures for composite structures. These procedures 
are naturally subject to revision as continuing studies enhance the understanding 
of composite failure mechanics. In conjunction with the theoretical develop- 
ments, there is the unsatisfied need for definition of suitable techniques for 
material property measurements. Accordingly, the program described here- 
in treated the above problem areas and the present report presents the results 
of investigations leading to: The definition of design criteria; the enhancement 
of the understanding of failure mechanisms’; the definition of unique design "con- 
cepts for composites; and the development of improved techniques for. com- 
posite property measurement. 

The section on "Design Criteria” describes the failure- criteria for lami- 
nates, which have been computerized for structural efficiency analysis. Fur- 
ther studies of the laminae failure mechanisms upon which this rs based are also 
described, including some treatment of time dependent behavior. The applica- 
tion of these criteria is treated in the "Design Concepts" section, wherein the 
utilization of the high uniaxial compressive strength of fibrous composites is 
emphasized. Columns, plates and panels of combined composite and metallic 
construction are designed and their potential is assessed. Concepts for im- 
proved material performance including three phase and isotropic three dimen- 
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sional materials are also defined. In the final section, ."Studies of Materials 
Properties Tests" the, 'results of the analysis of several current test tech- 
niques are presented along with suggested new methods for uniaxial strength 

measurements. 


' i* ..*c a.'-. . ;.n tbs tensl.;- load and -he re inf % ,c 

direction is incr s;r>3ci gradually from zero to ninety dcgim s, three 
primary failure modes can be expected to be encountered: first, zX small 
/.glee, lie failure of the filaments; seccud,* an Inte: -mediate angles. 
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DESIGN CRITERIA 


The parametric evaluations of fibrous composites for aerospace structures 

r' ‘ 

(e„ g. Refs.l and 2) have indicated the attractive potential of composite 

structures which are configured to achieve high stresses. When efficient 
stiffening arrangements permit the use of the high modulus composite materials 
high strength, the resulting structures are shown- to be substantially lighter 
than metallic structures. These studies emphasized the need for better 
definition of composite strength. The prediction of laminate strength from a 
knowledge of constituent properties is a complex undertaking and, in general, 
strength estimates most suitable fordesign are obtained experimentally. 
However, an analytical estimate is required' in the assessment of the potential 
of candidate composites which have not yet been brought to the practical fab-' 
rication stage. Indeed, the concept of analytically predicting composite pro- 
perties for use in a structural application analysis, is an essential part of the 
search for guidelines for the development of improved composite materials. 

The strengths of fibrous composites, perhaps to an even greater degree 
than their stiffnesses, are complex functions of the anisotropies associated 
with the uni- directional character of filaments. For filaments in one .(the 

loaded) direction only,, strength intension (Ref. 3), compression (e. g. Ref. 4). 

\ 

and shear, {e. g. Refs. 5, 6) have been related to. the properties of the consti- 
tuents. Further studyof these problems is described subsequently. Exten- 
sions to reinforcements in other directions have now been incorporated into 
the computer programs for structural properties' of anisotropic composite 
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elements (SPACE). The aim is to provide a strength assessment for prelimi- 
nary design studies. 

r ’ 

Numerous strength theories have been proposed; the validity of which can 
only be justified by experiments which are themselves complex and formidable. 
A comprehensive tabulation of strength theories has been presented in 
Reference 7. Among the theories proposed is the strength theory in References 
8, 9. In this theory, the strength of a fiber-reinforced composite, considered 
as quasi-homogeneous and anisotropic, is governed by a continuous failure 
surface of Hill (Ref. 10)’. Once the failure surface is determined, the strength 
of the composite body under any type of surface loading can be determined in 
a. straight-forward manner. This appears to be a reasonable approach for 
composites with elastic perfectly- plastic fibers and matrix'.' 

For most composites there is a vast difference between the strength of 
the filaments and the strength of the binder, and failures in the binder may 
be encountered as the maximum stress direction varies from the filament 
direction. In the simplest case, for example, of a unidirectionally reinforced 
composite, in tension, if the angle between the tensile load and the reinforce- 
ment direction is increased gradually from zero to ninety degrees, three 
primary failure modes can be expected to be encountered: first, at small • 
angles,, tensile failure of the filaments; second, at intermediate angles, 
shear failure in the binder; third, as the filaments b'ecome oriented mostly 
transverse to the load, tensile failure in the binder. These failure modes are 
essentially independent of one another. Particularly for the change from 
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tensile failure in the filaments to shear in the binder, there appears to be 
little reason to expect a gradual transition of the type that leads to a smooth, 

V 

"yield surface" for homogeneous materials. 

When the various failure modes are independent of each other for a com- 
posite having an oriented_ structure, the applicability of a continuous function, 
like Hill's anisotropic yield condition for a homogeneous material (Ref. 10), 
to represent a yield or strength criterion appears open to question. Accordingly, 
the approach utilized herein to determine strength criteria for composite 
laminates has been to 'determine separately the strengths for all possible 
failure modes. Thus, to a degree, a family of failure surfaces representative 
of the material will be generated, and the lowest of them for any loading 
condition will be the governing one. 

This approach is described in the following sub-section. This is followed 
by a presentation of the recent studies of the strength properties of a uniaxial 
laminate. Note that these latter quantities are required both for the dis- 
continuous and the continuous failure surface models. The analysis is developed 
first for composites where both fiber and matrix are elastic-plastic and a 
continuous failure surface is defined. These results are then specialized for 
rigid brittle fibers to generate strength values which can be used in a maximum 
sttess failure theory. 

Laminate Strength 

The strength analysis of a laminate of layers of uniaxial fibrous compo- 
sites utilizes the elastic analysi s, under given surface loadings, of the state 
of stress in each laminate layer considered as quasi-homogeneous, i. e. 
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locally heterogeneous^ but grossly homogeneous. If the surface loadings 
increase monotonically and proportionally, there will be a stage at which 
the stress in one (or more) layers of the laminate is at a failure point and 
the layer, being assumed to fail, is replaced by a new degraded layer having 
an assumed mode of degradation. As successive failure of constituting 
layers proceeds, a redistribution of stress among the laminae occurs and 
the slope of the load- deflection curve is discontinuous. The ultimate strength 
of the laminate is reached when all the constituting layers have failed. 

In each layer, the stresses referred to the principal axes of anistropy 
are computed. If the shear or transverse stress is equal to its corresponding 
yield stress, the lamina is considered to hold that stress level for those com- 
ponents and to have additional stiffness only in the fiber- direction. As 
successive failure ox constituting layers occurs, the entire load-deflection 
history can be traced until a failure in the fiber direction occurs or until all 
■layers yield, at which point the associated applied load is defined as the failure 
load. An illustration of the application of this approach is presented in Fig. 

1 where calculated 'stress- strain curves for two simple laminates of E-glass 
in epoxy are shown. For more general 'laminate configurations there will be a 

greater number of straight line sections in the stress-strain curve. 

/ \ ' 

For those applications where only the limiting stress levels - and not the 
entire stress-strain curve - are required, a simpler approach to the defini- 
tion of failure criteria appears reasonable. This approach is based upon the 
concept that the first 'departure from elastic behavior is a -most significant 



point on a composite material stress- strain curve. It is desired to keep 
actual stress levels below -this point in a fashion analag.ous to the use of the 
yield stress for metal construction. Similarly, at the ultimate stress level,, 
the transverse properties of the individual layers have generally deteriorated. 
Thus this level can be approximated by using a "netting" analysis with the 
uniaxial strength properties of the individual layers. The application of these 
principles in the definition of failure criteria is described below. -These 
methods have been incorporated as a subroutine of the SPACE computer 
program. 

The basic stress strain relations for the laminate are given by: 





where a ^ is the inverse of. the matrix a. 
For zero curvatures we have 


{.} =' La] _1 {n} 


“ d . - M -00. . 

<r. w = C.. s . 
i ij oj 

where & ^ are the stress components in the layer referred to the 
laminate axes. 

! /^\ 

C.. ' are the elastic moduli of the same layer referred to the same axes. 

(k) 

From these equations we may find the stresses, ex. , referred to the lamina 
principal axes from: 

^.(k) = ^ ,ok q. (k), no sum on ^ 

i ij 3 


where T is the transformation matrix for rotation of coordinate axes, 
k 

P is the angle from laminate principal axes to lamina principal axes. 

k 

To evaluate laminate strength,,, the stress components cr are evaluated for a 
laminate having its - laminae t thicknesses and the total load normalized. Thus, 
the total thickness is unity and the load vector for axial load, for example, is 

(1, 0, 0). With the stress components known, and the maximum strengths, 

(k) 

cr. x 7 defined for each of the stress components, (as in Refs. 3, 4, and 5) 

1U 

the ratios ^ (k) 

iu 

. (k) • 

< 7 . 

1U 
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are formed. The lowest ratio of this set of 3 n quantities (for an n layer 
laminate) is the material yield stress, cr . * 

The same analysis is now repeated for the case where the moduli E^ > 

G^» and are set equal’ to zero for all layers. Thus only E^, the modulus 

in the fiber direction is nonrzero and we have (for = v^E^J.a "netting" 

analysis. Here there are only n stress ratios and the lowest one is taken as 

the material ultimate str.es s, cr . 

u 

As an application, the yield stress of a symmetric biaxial laminate sub- 
jected to an axial tensile load was treated. The results for these laminates 
are presented in Fig. 2. .Each curve represents the results generated for one 
of the three' stress components and the lower envelope curve is the design' 
yield stress curve. 

Lamina Strength 

It has/been postulated (Refs. 8 and 9) that the strength of a unidirectional’ 
fibrous composite, considered as quasi-homo^eneous and anisotropic, can be 
represented by a surface having the form of Hill's Generalized von Mises' 
Yield Condition (Ref. 4), namely: , 


Zi ( T ij) F ( T ZZ~ r 33) +G ( T 33 _r il) + H ( r il” T 22.) + 2LT 2 


2 

23 


■ (1) 


+ 2M T Z + 2NT' 2 ' = ' 1 

31 12- 


where the coefficients F, G, H, L, M, N are parameters characterizing the 
state of anisotropy and T are components of the stress tensor referred to the 
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principal axes of anisotropy x^ x^ x^ where the x^ axis is parallel to the 
fibers. * In general, the yield condition ( 1) can be represented by a surface in 
a stress space. A stress point within the yield surface represents a "safe" 
stress state. Yielding can occur if the stress point is on the yield surface. 
Since unidirectional fibrous composite layers in laminates are thin compared 
with their lateral dimensions, only and T ^ are considered non- zero . 

if the x^ axis is along the thickness of the layer. Furthermore, since uni- 
directional fibers are randomly located in a composite layer, it can be assumed 
to be transversely isotropic. Then the yield condition (1) reduces to 



where X^ , and X^ are the normal yield stress in the direction of the 
fibers, the normal yield stress in the direction transverse to the fibers and 
the yield stress in axial shear of the composite, respectively. These are the 
three basic strength characteristics of the unidirectional fibrous composites. 
Once these are known, the yield condition (2) can be employed to determine 
whether a combined state of (plane) stress can cause failure of the composite. 

In what follows, effort is made to evaluate the quantities, X^ , X^ and 
X^ analytically in terms of the strength and geometry of the constituents. 


# Henceforth, unless otherwise specified, i, j = 1, 2,- 3; Summation on 
repeated indices is implied. 
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Limit Analysis of Unidirectional Fibrous Composites 

The composite material under consideration consists of a 

■(' 

relatively soft matrix material in which stiffer fibers are embedded. Initially 
both materials ar.e assumed to be elastic-perfectly plastic and satisfy the. von, 
Mises' yield criterion. As shown in Figure 1, referred to an orthogonal 
Cartesian co-ordinate system, a typical unidirectional fibrous composite is 
taken to be a cylinder with rectangular cross-section. Circular fibers run- 
ning from base to base of the specimen are in x^ -direction. Limit analyses 
of such a specimen with various arrangements of both elastic brittle and elastic- 
perfectly plastic fibers in an elastic-plastic matrix under various types of 
surface loading will be described. 

Elastic -plastic Constituents 

In this study only the "random .array" geometry is considered:, circular 

fibers of various diameters are randomly located in the specimen. Each of 

them can be surrounded entirely by a concentric cylinder of matrix material. 

A cylinder consisting of a fiber of radius r^ and the outer matrix:- shell of 

radius r^ is called a composite cylinder. It is assumed that a constant 
r 

£ _ f can be chosen so that the composite cylinders are non-overlapping. 

‘ \ 

The entire specimen can then he considered as an assemblage of composite 
' cylinders and the ramining matrix volume. The lateral boundary of the spe- 
cimen may touch' or cut through some fibers. In both cases the associated 
composite cylinders are "incomplete". Since in practice fiber diameters are 
very small compared with the transverse dimensions of the specimen, the 
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total number of fibers in the interior of the specimen is much larger than the 
total number of those possibly on the lateral boundary. Hence, the total 
volume of "incomplete" composite cylinders is much smaller than that of the 
"complete" cylinders. 

If V, V ^ and denote, respectively, the total volumes of the specimen, 
the composite cylinders and the remaining matrix in the specimen, the 
following obvious relation holds; 

V = V 1 + V 2 (3) 


In the case where the entire specimen is occupied by composite cylinders, 
V 2 = 0. Following Hashin and Rosen (Ref. 11) this distribution is called the 
"random array". 

The von Mises" yield criterion which the fiber and matrix materials are 
assumed to. obey has the following form (Ref. 12); 


S.. 

U 

2 


S.. 
_J 0 


(4) 


where S. . are components of the stress deviator and k is the yield stress in 
simple shear for the fiber material (denoted by k^) or for the matrix material 
(denoted by k^). 

Under the conditions of plane strain perpendicular to the x^-axis, von 
Mises' yield criterion (4) reduces to 


T - T 
22 33 


2 . 2 
+ 4 T 
■ 23 


4 k 


(5) 
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where , T ant ^^3 are components of the stress tensor in the transverse 
plane . 

In order to evaluate the three basic strength characteristics (Refs. 8, 9) 
of the composite specimen shown in Figure 3, the specimen is subjected to 
axial shear stress transverse tensile stress and longitudinal tensile 

stress , respectively. The upper and lower bound theorems of limit 
analysis of plasticity (Refs. 13, 14) will be used to obtain bounds for the limit 
loads T ^ , T and T ^ ^ w hich represents the lamina strengths'. 

Case 1 . Axial shear stresses T applied on the boundary of the composite, 
specimen. 

According to the lower bound theorem, a uniform shear stress field 
can be chosen as the statically admissible stress field. Since can nowhere 
violate the yield condition (4) for both matrix and fibers, it follows that a . 
lower bound for the limit load T ^ is 



For upper bound construction, a kinematically admissible velocity field 
is chosen as follows: 

(a) In (and thus also on the boundary of the composite cylinders) and 
on the entire lateral boundary of the specimen. 



where is a positive number. 
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(b) In any composite cylinder, referred to a local cylindrical 
' polar coordinate system, * 


- Yj. r cos B 
Yj Xj cos 0 

- Yj sin 0 


for 0 r s 


V ^ 

~ C 1 - r cos 

1 -P r 


cos 0- 


for r . 22 r s; r. 
f b 


- Y j Xj sin 0 


* The velocity field in any "incomplete 11 composite cylinder is defined by 
solving similar elastic displacement boundary value problems for the 
"incomplete" composite cylinder. However, since the volume of the "in- 
complete" composite cylinder is small, the difference between their actual 
contribution to the dissipation function and that obtained by treating all • 
cylinders as "complete" is negligible. This approximation is implied in 
the subsequent analysis wherever a similar situation arises. 
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The velocity field (7) in a composite cylinder is the elastic displacement 
solution to the bo'undary value- problem with the boundary condition (6) pre- 
scribed. The problem is the same as that which was formulated by Hashin 
and Rosen (Ref. 11) with the modification that for 0 <r ^r^ , the velocity field 
is associated with rigid body motion. 

'With^this velocity field constructed for the entire specimen, the dissipa- 
tion density function and the rate of external work done can be obtained to 
yield an upper bound for T ^ (Ref. 5). 


2 IT 


llUJL = f f 


P 



1+ 2 cos 0 d 8 dR (8) 
R 


The above expression is for "random array" in which = 0 and the fiber 

^ 2 
volume fraction v^ = P • 

The integral in (7) is "integrated numerically for different fiber-volume 


fractions. The result is shown in Figure 4 where 




v in (8) is plotted 


V 


as a function of v^O*, v f «= 1, v^ = P ). Note that particularly in Figure 4, 



From the above result, it is concluded that is at most about 27% above 
and at least the same as the yield stress in shear for the matrix. 
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Case 2 „ Transverse tensile stresses 7 ^ applied on the boundary of the 
composite specimen. 

For lower bound construction, a uniform tensile stress 7-^ throughout the 
specimen is chosen as a statically admissible stress field. Since the von 
Mises' yield condition (4) can nowhere be violated in the specimen, the lower 
bound associated with the constructed statically admissible stress field is 

• {’£ ) ■ 151 

For upper bound construction, a kinematically admissible velocity field is 
chosen as follows: 

(a) In V 2 (and thus also on the boundary of the composite cylinders) and 
on the boundary of the specimen. 



where is any- positive number. 
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(b) In any composite cylinder, referred to the local coordinate system, 



x? + 3 x 2 


•2 T J *3 + r { <‘ x 2 + 3X 3 5 + A < K ' < > 

•5 , A ■■■' 


2.2 2 

*3 


2r„ 


2 r 


3X 3 + X 3 + A, !>ili>, + A^f <‘ x 2 +x 3> 


J 


2 r J J 


for r. 


( 11 ) 


0 

V 2 2 

J2 x B + B, X 2 + 3x 3 
2 2 1 2 . 




2r, 


- y_ 3x 2 + x 2 

J2_x B + B, 2 + x 3 
2 o i 2 ^ 


for 0 < r r 


2r 


where 


A i=^ 4 

a 2= i 


: (l-Ti)(l+ri) g 4 - 3(l-ri)(l+r)) ^ + 


2 i 


„ 2 . „2 


1} 


4 (i-rD(i+n)(i-in g 

A 3 = ^ ' |^( 1-r l) '[ ( 1 — rj) P 2 + M 

a 4 = h { 4(1 - T '> [-211+ (1+n) (1+ )] 

B ! = £ jp (1-Ip P 4 - 6 (1-T)) p 2 + 2(1+11 ) l 


B 2= C 


e (i-n) g 2 (i-p 2 ) 


and C = (l-n') 2 g 6 + 4(l-ri)(l + n) g 4 . 6 (I-t, ) (1+n ) g 2 + 4 (1-^) (1+ n) + 


(l+T)) 
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The velocity field (11) in a composite cylinder is the incompressible 
plastic displacement solution to the boundary value problem with the boundary 
condition (10) prescribed. The problem was formulated by Hashln and Rosen 
(Ref. 11) in evaluating the effective plane strain shear modulus of elastic 
fiber- reinforced composites. The constant rjin (11) considered as the displace- 
ment solution to the elastic problem is the ratio of the elastic shear moduli 
of'.the fiber and matrix materials. However, it can be considered as merely 
a parameter when (11) is used as a kinematically admissible velocity. 

For the case of the "random array", an application of the upper theorem 
gives an upper bound . • 


( T 22 )u X + I .(12) 

2k b \ 


where 




4 2' 2 < 1+T 1) ' 2 2 

8(1- Tl) p 4 -6 (l-ri) r + p 2 + 12(l-n) (1-P ) R- 

c 


dR 


and 


I 


2 


1 1 2ff R 

1 r j f c" 7^ (Tl.,'p, R) + i 2 (n, P,R) cos e - d0 dR 


in which 

h 


#,(n, P,R) =. U (l-ri)(l+ri) : p 4 -3 ( 1 -T|) (i+n) p 2 + 6 ( 1 -T 1 )( 1 +T 1 )( 1 -P 2 )R 2 


(1-^)¥.|- 3 


.R 


R L P 


(I+Hi 

o2 


+ (i-ti) r 


+ 2f 


-2n+('i+n)(.i+ 

bO 
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® z Cn,.P,R)= 2 ;4(l-Ti)(l+n> P -3(1-Tn(l+Ti)p 2 + i^^+6(l-r|)(l+Ti){l-p 2 )R‘ : 


(x -rp P .3 


— { i --n) p 2 


+ 2 p 


-2T] + (1+T1)(1+ - 6 ) 


and P = v^. 


For any finite k, , and for r| -* .» , (12) becomes 




1 2ff 


2 1 ^ ff(l-P 2 ) 3 


J J B. (P, R) +''1' 2 (P,R) cos B d0dR (13) 


$ o 


where 


4 f 2 

V P * & )= -V O 4 +p 2 +l)- - I- (f3 2 + 1) p_ 
R L • 2 R 2 


-v 2 


W - 4 ^ 2 ± L ) + 3P 2 r 2 


n 2 


and 


^(P.R) = 

2 R 2 


i 4 + P 2 + l)- |(P 2 +1) *— 

2 R 2 J 


(4g 4 + P 2 + 1) 3 p 2 R 2 

2 


The right hand side of expression (13) was obtained in Ref, 5 'for 
an upper bound of the limit load (normalized with respect to k^)-for transverse 
shear stresses applied on composites reinforced with rigid fibers. In Figure 

(r h ) - 

5, \ 22 / v is plotted as a function of v according to (13). It is a monotonically 

2 k^ * • 

increasing function of v. with lim fr,,') TT \ , , lim /r „ 

f v -> o \ILLE = 1 and [ZZ ) U = » . 

f 2 \ ' f ' 2 \ 

On the other hand, for any finite T), y can be plotted as a function of 

2k b 

by evaluating numerically I, and I, in (12). For example, for r)= 100, T 22 'u 
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is also a monotonically increasing function of v 


HV° 


22 j U 1 x and 
2k b 


v 1 \rr - — = t — ) *» which is shown in -Figure 3, with f = 10 . It is 
' f k b , / - Jc . 

r ■ D 

observed that for fiber volume fractions smaller than about 75%, the upper 
bound for the limit load i-s higher that for T[~* » but-for higher fiber volume 
fraction, the opposite is true. Therefore, for .any fiber volume fraction between 
0 and 1, the lower value of the bounds obtained from T|=: 100 and T)-*« will give a 


better estimate of T' and it is shown in Figure 5 in solid line as a 


function of 


It is interesting to note that for T) = 1, (12) reduces to the following 
simple form 


tsL)= l + v (-L - i) (14) 

2k b . f W / 

which is commonly known as the "rule of mixtures". From Figure 3, it is 

seen that the straight line represented by (14) will be higher than the chosen 

(* h \ . - 
curve for v 22 / U for 0 < 1 . Therefore, it is concluded that the "rule of 

2 \ . ■ ; 
mixtures to be used in' this case would overestimate the composite strength. 

In Figure 5 , ( T 22 } L» obtained from (9) is also shown. The difference 
■ . •. ' 2k b 

between the upper and lower bounds for v^-'O is due to the fact that in the 
upper bound construction, a plane strain velocity field is used as' a kinemati- 
cally admissible velocity field. 

Case 3. Longitudinal tensile stresses Tj applied on the boundary of the 


composite specimen. 
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The study of longitudinal strength of unidirectional fibrous composites is 
extensive. Various models and failure mechanisms have been proposed in 
the literature. Kere p the yield strength in axial tension is obtained by, the 
construction of very simple velocity and stress fields, for composites having 
elastic-plastic fibers of uniform strength. 

For upper bound construction, in the entire region of the composite spe- 
cimen, a kinematically admissible velocity field is chosen as 



An application of the upper bound theorem gives 

( r n) = T f v f + Vv ■ 



■where T and T are the yield stresses in tension forv the fiber and matrix. 


respectively and v = 1 - v f . - 

Relation (15)- is known as the "rule of mixtures" for the predication of the 

7 T L \ 

strength of a composite. If T and T are given, V 11 / U is a linear function of 

‘ ' T b • 

fiber-volume fraction v^ . 
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For lower bound construction, if the applied stress T ^ is assumed to be 
uniformly distributed on the boundary of the composite specimen, then a 
lower bound ^ is equal to . However, if we assume 'that the tensile 

stress is not uniformly distributed on the boundary surface and we are 
only interested in the average stress intensity on the boundary surface that 
causes failure, then a higher lower bound can be obtained. In fact, a statically 
admissible stress field can be chosen as follows: 

In the region occupied by fibers 
T 0 0 

O' 0 ' 0 

0 0 0 

In the region occupied by the. matrix, 

0 0 
0 0 
0 0 

It is obvious that the stress field expressed in (16) and (17) satisfy the re- 
quirements to be statically admissible. 1 

The average traction corresponding to this stress field is therefore 
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this type of composite, {i. e. elastic-plastic fibers and- matrix of uniform 
yield strength). 

From the above results, it is observed that obunds for T ^ are far apart, 
especially for high volume fractions. Further effort should be made to im- 
prove the bounds in order to have a better estimate of the limit load. For lower 
bound construction, uniform stress distribution used as statically admissible 
stress fields can only give lower bounds which are the corresponding matrix 
yield stresses. In order to obtain higher lower bounds, one has to assume 
applied tractions to be non- uniformly distributed on the boundary surface 
according to the properties of the fibers and matrix. Equilibrium stress fields 
can then be constructed in equilibrium with the applied tractions. Then the 
lower bound theorem can be applied to obtain higher lower bounds. To decide 
the distribution of the applied tractions on the boundary and to construct, an 
equilibrium stress field in the body is not at all easy in general. The 'success 
in the construction of such a statically admissible stress field for Case 3. is 
due to the simplicity of geometry and loading conditions. 

Brittle Fibers 

In the previous section, the strength of unidirectional fibrous composites 
has been evaluated by obtaining bounds for the basic strength characteristics. 
The theory assumes the existence of a continuous failure surface. This hypo- 
thesis appears reasonable for composites with elastic-perfectly plastic 
fibers and matrix such as metal fibers. and matrix. However, for contemporary 
high strength fibers embedded in epoxy resin matrix, this "elastic-perfectly 
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plastic model” may seem inadequate. Instead, it seems more suitable to 
assume elastic-brittle fibers and an elastic-perfectly plastic matrix which 

•f' ' 

obeys the von Mises' yield criterion. Upper and lower bounds have been ob- 
tained for the following types of surface loadings: 

1. Tractions equivalent to a uniform shear stress T ^ applied on the en- 
tire'boundary surface of the specimen. 

2. . Tractions equivalent to a uniform shear stress applied on the 

entire boundary surface under the conditions of plane strain. 

3. ' Tractions equivalent to uniform uniaxial tension applied on the 

entire boundary surface under the conditions of plane strain. 

4. Tractions equivalent to biaxial uniform tension and T 

Ac 3 3 

(r 22 4 r 33 ) applied on the entire boundary surface under the conditions 
of plane strain. 

5. Tractions equivalent to combined in-plane shear and transverse 
tension (T^ j = a ) applied on the entire boundary surface. 

Since the method of analysis is similar to what has been presented in the 
preceeding pages in obtaining bounds for the basic strength characteristics, the 
details of analysis will be presented in Appendix A. However, results are 
summarized as follows: 
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for the "random array" geometry. • The above expression for \ 12 /U is the 

' ' h 2 

same as- (8) which is shown in Figure 4 as a function of v^ (0 < ^ ^1, = P ) 


since the same velocity is used here. 

For arbitrary geometry of arrangement of fibers 

which is higher than ( T 12 ) U given by (8), 

I L\ 


\ . 


1 t-'U _________ 

U _ - f f Vf (p 7 r) + VP.R) cos 8 d 6 dR 

- ird-pV I I 1 


•where 


yP.R) 


£ fs 2 


( 4 >W 


±JL)+ 3 P 2 R 2 


and yp, r)=^y p 4 +p 2 + i -|(p 2 + 1) -C- 


+ 3 p 2 R 2 


for the "random array" geometry. The above expression for ( T 23 ) U is the 

i v • i.i at ~c j. l ' „"V } ; j “he uppc-r do.. / id cypres , o,- 7 

same as (13) which is plptt e 4 . in^Figure ..5 , as ^ipnction.o| sirer-^n of rho 

idu mvjoiiad. Omls ,s dae to ihe fa'cL thou c^ody motion m 

11111 ( T 23 )u _ 1 , 1 ™' \ 23 /U 

-,ys Vt:° -hri^TStS Le idnemalSM'y Vt^rr =H3;ce, r.gwt. 

since the velocity field chosen here can be obtained from (10) and (11) through 
an orthogonal transformation together with the limiting process, of T)-* ® 
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and 


for 


and 



4. ■ 

( T 22)u 

(rf)u = 2 V 


,( T 22 L )u - 1 

*-3 CO 

CM 

11 

^ CO ‘ 

5.. 

II 

»-3 CM 

1 . 



/ 1 2 
Vs + “ 


tllli = 

O' 




a a O'. 



For the 

case o£ the 

"random array" \ 


( T 22.)u _ 

I, 

min 3 


■s 

1 

< > J + OdCO 



a I, 

xniiji 3 


> 

■ M i 

1 ' ■ 2 + OiW. 


( 19 ) 


( 20 ) 
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where 


I 


3 


7T(1-P 2 ) 3 *' 


1 27 r 



H (to, P,R)+0 2 (a),P,R) cos B+jyP.R) cos 2 6 


n i (w, P, k) =4 


i 4 + p 3 + 1)-|(P 2 +1)-|(#+1)^1 
2 2 r_ 


3 A 2 - <*£±£± 1 ) 


-,2 


+ w 2 (1-p 2 ) 4 (1 +■ 


a (W, P, R) = 2 w 2 (i-p 2 ) 4 Jb 
2 R 2 


and (P, R) = 8 • 


(p 4 + p 2 + 1) - 1 (p 2 + 1) 


3 p 2 R 2 . (4 e 4 +• £ + 1) 


Numerical calculation is performed to obtain and (V from 

\ 22 /jkJ \ 12 JX} 

(20) for different values of P and Ct . The results are summarized in Figure 6 
2 

in which {3 =0.8 is the highest fiber volume fraction shown. The dotted line 

represents (19) which gives the lower bound for ^^for any (3 . 

It is worth mentioning that for this model of elastic -brittle fibers embedded 
in elastic -perfectly plastic matrix, in the upper bound expressions, only the 
matrix strength and fiber -volume fraction appear — the brittle strength of the 
fibers is not involved. This is due to the fact that rigid body motion of fibers is 
always incorporated into the kinematically admissible fields-. Hence, regions 
occupied by fibers in the composite body contribute nothing to the dissipation 
function. The result is reasonable since in reality, elastic moduli and brittle 
strength of fibers are much higher than the corresponding elastic moduli and 
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strength of the matrix material. On the other hand, the longitudinal strength 
along fiber -direction of fibrous composites under tension or compression cannot 

V* 

be obtained by limit analysis of plasticity because of the presence of the high 
modulus, high strength elastic -brittle fibers. Instead, a statistical failure 
theory was established for tensile strength [ 3 ] and a fiber -buckling model was 
postulated for compressive strength of fibrous composites [4] and [l5], These failure 
mechanisms are based on experimental observation. Furthermore, as men-' 
tiohed previously, for this mathematical model of composite material, the 
existence of a smooth yield surface that governs the strength of the material is 
still an open question. Based on experimental investigation, the failure mech- 
anisms for individual cases seem to be different and independent of one another. 
Accordingly, the approach to the determination of the strength criterion for uni- 
directional fiber -reinforced composites as well as fibrous composite laminates 
is as described earlier. Instead of using a continuous yield surface, all possible 
failure modes are considered. To determine which one is dominant for a parti- 
cular type of loading condition can sometimes become a tedious task especially 
for fibrous laminates under various types of design loading. However, this 
approach has the advantage that the entire load-deflection history can be traced 
until complete failure occurs. With modern high speed digital computers avail- 
able, a systematic strength analysis can be programmed to obtain accurate 
results. 

For certain composites, it is possible to utilize, as an alternate approach, 
the application of limit analysis methods to the laminate as a whole. This is 
treated in the following section. 
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Application to Laminates 

It is interesting to note that the kinematically admissible velocity fields 
constructed to evaluate the upper bound for the three basic strength characteris- 
tics can be used to obtain upper bounds of limit loads for in -plane loading applied 
on laminates . 

To demonstrate the method, consider a laminate composed of uni- 
directional fibrous laminae subjected to in-plane shear .stress * ^^e 

L 

.problem is to find an upper bound for the limit load T . A typical constituting 

12 

th - ( k > 

layer, the k layer, in which fibers are all running in x^ direction, is shown 

in Figure 7. Depending on the orientation of fibers in the layer, the principal 

axes x x can be defined by an angle 0 ^) measured from the lami- 

12 3 

nate axes x x^x . 

• 12 3 

For upper bound construction, a kinematically admissible velocity .field 

is constructed in the following manner: 

The velocity field (P , p , P ) = (0, y, x, , 0) is assigned in V- 

1.2 3 c. ~ ' 

of each layer, referred to x^ x^ x^ axes. In the k^ 1 layer referred to its prin- 
« cipal axes x^) x x^ > the above velocity field is transformed into the 
following form: » 


Y lS ine< k > cos e (k) Xi < k )= Yi sin 2 0< k > X2 «' 


Yj, cos 


2 < k) Xl (k) -Y 1 sine( k )cos0 (k) x 2 < k > 


The right hand side of equation (21) can be decomposed into four parts 
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(22) 


*S {k) = + + + 


1 * 2 


3 


where 


M M. 
.1 


Y ^ cos 


2 0 (k) x (k) 


- (k) 
^2 


“V c . 0 (k) a (k) (k) 

— sin o cos t) x^. 


Y 1 . fi (k) fi (k) (k) 
— sin 9 cos 0 x 3 


Y 1 sin6 (k) cofl 0* k) Xj (k) 


-(k) 

U 3 


1 . ft (k) ' fi (k) (k)‘ 

sin 0 ■ cos 0' 


-Y 


1 . fl (k) 0 (k) . (k) 

2 ~ sm 0 cos 0 x 3 


and 


- (k) 

V 


-Y, sin 2 9 (k) x (k) 


Except for multiplicate factors and reference coordinate system, . 
u.^ \ u , (k >, and u,^ are the same as those velocity fields constructed 
in for the three different cases of loading in the proceding section. 
Furthermore. 


30 - 



(k) = 


, . 2 0 (k)- (k)' 

i-Yj sin 6 W x 

0 

■?' 

0 


. . 2 .(k) (k) 

= <-Yj sm 8' _x ' 


V • 2 fl( k > v (k) 

-Y^ sm o 

v • 2 fl< k > < k > 

sin o x^ 


(23) 


— (k) 

where the first part of the decomposition has the same form as u. ' while the 

second part represents a rigid body motion of as a whole. 

The velocity field chosen in is' also composed of four parts:' The first 
— (k) 

part is due to u , which is the same as (7.) with the modification that Y, 
cos^ 6^. The second part is due to u ^ -which is the same as (11) with the 

/ -4^ 

(k) (k) 

modification that Y^ is now replaced. by -Y^ sin 0' 7 cos 9^ . The third part 

— /jk) (k) ■ ’ 

is due to u_ which is exactly the same as u ' in Y . The fourth part is due 

— (k) 

to u x 7 which can be obtained in the same way as in Case 1 in the preceding sec- 
-> 4 

tion. However, a moment's reflection reveals that the strain rate associated . 


with u 


r fr> „ 


2 4 


n (k) 


Y^ cos 8 V 7 is replaced by -Yj sin 6 V 7 since the additional rigid 
body motion contributes nothing to the value of the strain rate. Now, after 
the kinematically admissible velocity field is constructed as described above, 
the associated strain rate and dissipation density function can be calculated 
without .difficulty. Moreover, the rate of external work done can also be cal- 
culated so that an upper bound for t ^ can be obtained. * 

As an example, consider a laminate composed of two unidirectional 
fibrous composite layers of equal thickness with 0^ = 0 and 9^ ^ . 

In this simple geometry of lamination, the kinematically admissible 
velocity field constructed above reduces to a very simple form. 

For layer 1, the velocity field (21) reduces to 
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( 24 ) 


For layer 2, the velocity field (21) reduces to' 



(25)' 


Then, following the principles described above, it is easy to obtain the velocity 
fields and the associated- strain rate fields for' both layers. After some mani- 
pulation, it turns out that 


(T L ) 

' T 12 'U i ■ _ 1 

k b tr(l-p 2 ) 


1 2tt 

I / 


R 



cos 0 d9 dR 


for the "random -array" geometry. 

The above equation is exactly (8), which means that the upper bound 
* s not higher for this type of cross-ply laminates than for the uni- 
directional fibrous composite. 

. Similarly, it is easy to show that the three, types of velocity f-ields con- 
structed in the preceding section can be used to construct upper bounds for 
limit loads for laminates subjected to any in-plane stresses. 
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Tensile Strength 

The high-strength, high modulus fibers which are of interest for 

j •. 

use in composite materials are generally brittle, having tensile strengths 
that must be characterized statistically. Any theory for the tensile strength 
of composites containing such fibers must take into account the dispersion in 

their failure stress levels in order to have any relevance. 

(16) 

Parratt noted the dispersion in fiber strength and suggested that 
failure of a fibrous composite subjected to tensile load occurs when the 
fibers have broken up into lengths so short that any increase in applied load 
cannot be transmitted to the fibers because the limit of interface or matrix 
shear has been reached. 

A theory has teen presented (3) for the failure stress of composites 
containing continuous, uniaxially- oriented, brittle fibers in a ductile matrix 
loaded parallel to the fiber direction. This theory predicts that due to the 
distribution of flaws or imperfections in the fibers there occurs a series of 
randomly -distributed fiber fractures as the applied stress level is increased. 

It is argued that a portion of the broken fiber in the vicinity of the 
fracture is ineffective in resisting the applied load. Assuming that the stress 
in a broken fiber is uniformly distributed among the unbroken fibers in the' 
cross-section and that this over stress acts over a length equal to the- "in- 
effective" length, it is predicted that failure occurs when a weakened' cross- 
section cannot sustain an increase in load. In effect, the theory pre- 
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diets that the composite has a strength equal to that of a bundle of fibers 
whose length is the "ineffective" length. 

The present study considers fibers having a statistical strength dis- 
tribution resulting in fractures at various stress levels as the'applied load 
increases.. It is assumed that the stress in the broken fiber is distributed 
unevenly to the other fibers in the cross-section which has a length equal to 
the ineffective length. As a first approximation the effect of this overstress 
is presumed to affect only those fibers adjacent to a break. Failure is as- 
sumed to occur due to an increasing probability of fracture in the fibers ad- 
jacent to a prior break. 

Description of the Model 

The model consists of a two-dimensional composite of length L 
consisting of a ductile matrix in which are imbedded N continuous brittle 
fibers whose orientation is parallel to the applied tensile load. The fibers 
are considered to be composed of layers of length 6 . The total number of 
layers being M = L / 6. (See Figure 8, ) 

• The quantity 6 represents some length’ over which the stress is per- 
turbed in the area of a fracture. It is variously referred to as the ineffec- 
tive length or twice the transfer length and several formulae have been. proposed 
for its evaluation. Two of the definitions(3), (17)are based on an elastic shear - 
lag type of analysis* The axial stress in a broken fiber is found to be 
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( 1 ) 


a^(x) = CTq [l +; sinh r)x - cosh r|x] 


where 




1-v 


1/2 



b 

r 

f 

V f 

X 


•fiber end. 


= Young’s modulus of fiber 
= Shear modulus of binder (matrix) 

= Fiber radius 

= Volume fraction of fibers 

= Dist ranee from end of broken fiber 

= Extensional stress in the fiber at a large^distance from the. 


It should be noted that in Reference 3 .the factor of 2 in the expres- 
sion for r\ was incorrectly omitted. The first author -defines 6 ;as the value 
of x for which the stress in the fiber has reached 90% of O' » the stress at a 
long distance from the fiber break. On the other hand, Friedman defines the 
ineffective length by means of an approximate step -function stress distri- 
bution which has the same average stress as the distribution of Equation 1. 
Furthermore, this author includes the ineffective portion on both sides of 
the break whereas the first definition uses only one side. 

The expressions for the two ineffective lengths discussed above 
are, respectively 
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If the shear stress between the matrix and fiber ‘is assumed to be 
constant, as in the case of plastic flow or frictional shear stress due to inter- 
face failure, the fiber tensile stress becomes constant at a finite distance 
from the fiber end. This distance is called the transfer length ^ by some 
authors and is given by the expression 



where T is the (constant) plastic shear stress of the matrix, or, if there is. .• 
interface failure, the frictional force between the matrix and fiber." 

The ineffective length can also be determined experimentally by a 
photoelastic examination of the stress patterns in the vicinity of a fiber break. 


In thi£ investigation’ the experimentally determined ineffective length is used 
. • ’ \ 
where available. • 


Statistical Analysis of the Model 


} • 


It is assumed that the strength of the population, of fiber elements of 


length 6 can be characterized by a cumulative distribution function F(cr). That 
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is, the probability that an arbitrary element has a failure stress level less 
than, or equal to 0 is F(cf>). 

The expected number of elements in the composite that will fail under 
fiber stress a is 

E = MNF(a) ..., • ' (4) 

This expression includes the possibility of further fractures of adjacent 
fibers. If the composite does not fail, the stress in the broken fibers is dis- 
tributed to the other fibers in the cross-section in a complex manner. Asa 

(181 

first approximation, Hedgepeth^ ^used a shear lag analysis to determine the 
average stresses in fibers adjacent to an arbitrary number of broken fibers. 

He considers an infinite two-dimensional array of fibers subjected to 
tensile load parallel to the fiber direction which is uniform at a great dis- 
tance from the flracture area. The ratio of stress in the two fibers adjacent 
to a run of r broken fibers to the uniform applied stress at infinity is, for a 
static stress distribution 


_ 4-6-8 *(2r+2) 

r " 3*5*7* • (2r+l ) 


(5) 


Hedgepeth calls a stress-concentration factor, but in this paper it will be 
referred to as an overstress -factor so that it will not be confused with stress 
concentration factors found by an ,r exact M analysis of the stress distribution. 
For the case where r fibers break simultaneously the author demon- 
strated that the ratio of the maximum dynamic stress to the static stress in 
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fthe fibers adjacent to the break increases from 1.15 for r = 1 to a limit of 

i- 

1-. 27. The values for .,r = 2'‘and r = 3 are 1. 19 and 1. 20 respectively. 

• In the present paper it is assumed that the overstress in two fibers ad- 

jacent to the broken ones exists over the entire ineffective length. Therefore, 
the probability that an element adjacent to r broken elements, will fail is, 
approximately, the probability that its strength lies between a and K a . 

This probability is equal to 

. F(K a) - F(CT) (6) 

This approximation is justified, as will be shown later, by the fact that com- ■ 

posite failure occurs for small values of F(c) and therefore the probability of 

having adjacent fractures because both fibers have strength-less than a, which 

2 

probability is proportional to F(a) , and the probability of interaction of frac- 
ture groups is relatively small. 

Given that a single element is broken the probability that one of the tw.o 
adjacent fibers will break is 

2 . 

p 2/ l = 2[F (K iCT) - F(a)J - 2 [F (K i4 - F{a)] . (7) 

The. probability that both adjacent fibers will' break simultaneously is 

P 3 /l= [F ( V) -F(a ) ] 2 (8) 

.It should be noted that each of these expressions does not exclude the possi- 
bility of further fractures. 
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It is now assumed that only the two fibers immediately adjacent to a break 
are subjected to an overstress and that all of the remaining fibers in the cross- 
section have a stress level equal to the average stress cr. If one of the fibers 
adjacent to a single fracture breaks the fibers adjacent to the two broken 
fibers are subjected to a stress level K^. One of these overstressed fibers 
was previously exposed to a stress level K^o while the other saw only the 
average stress CT . The probability that one of the two fibers will break is 

?3/ 2 = [F<V)-F(Ka)] + [F(K 2 a,-F(a)] (9) 

-2 [F(K c) - F(K 0)] [f'(K 2 <J) - F{a)J . 

The probability that both fibers will break simultaneously is 

P 4/2 = [F(K 2 CT) - F(Ka)] [F(K 2 a) _ F(a)] (10) 

If both fibers adjacent to an initial fracture break there will be three 

broken fibers in a row and therefore the two fibers adjacent to this group, 

which were previously at a stress level a, will be subjected to a stress K^a . 

Again it is possible for -one or two of these fibers to break, and so on. 

By now the process and complexity of thfs problem should be evident so 

that the expressions for further fracture probabilities are presented without 

\ 

discussion. Note that that there are two paths by which a state of three 
broken fibers can be reached from a single break; A) by the simultaneous 
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breaking of both fibers adjacent to the initial break, B) by the successive 
breaking of two fibers. 1 

The expression p. .. repre'sents the probability of having i fibers 
1/ J z 

broken given that j are already broken. The letter z represents the particular 
path if there is more than one. The letters A and B refer to the paths des- 
cribed above. 


P 4/3A' = 2 [ F(K 3 ff) ‘ F(G) ] - 2 [ F < K 3 a > - F < t7 0 2 

P 5/3A = [ f <K 3 cj) - F(a).] 2 (lla-d) 

P 4/3B = [ F(K 3 ff) ' F( V>] + [ F(K 3 a) ' - F(CT) ] - 
2 [f(K 3 ct) - F(K 2 a)J [f(K 3 ct) - F(a)] 

P 5/3B = [ F < K 3 ff) - F(K 2 a)] [f^o) - F(a)]. 


Probability of Cumulative Fractures 
Each of the E^ fracture sites distributed throughout the composite is a 
nucleus .for further fiber breaks because of the overstress in adjacent fibers. 
The probability that an element will break followed by the fracture of at least 
one adjacent element is 

P 2 = F(a) (P 2/1 + P 3/1 ). <«) 

Therefore, for the composite as a whole the expected number of groups of 
two or more -broken fibers is 
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' (13) 


E = — — — , 2 
2 2 


and the associated probability of having at least one such group is 


P 2 = H1-}p 2 ) M(N - 2 > 


(14) 


■where the factors of 1/2 are introduced to account for the independence of or- 
der of two fractures. 

The probability that an element will fracture followed by the breaking of 
at least two other fibers in a row is equal to the probability of at least two 
fractures less the probability that two will break without further fractures 
occurring. This probability is 


P 2" F(t7)P 2/l (1 - p 3/2- p 4/2 ) 


(15) 


The probability of having at least one group with three or more frac- 
tures is 


1 M(N-2) 

p 3 -i- ( i- T p 3 ) 


( 16 ) 


and the expected number of such events is 
E 3=T P 3 M(N - 2) 

The analogous expressions for groups of four v or more fractures is 


P 4 = P 3 -F(a) [p 3/1 ( 1 -P 54/3A ) + P E/l P 3 /2 (1 " P 54/3B ) ] 


P 4 


= 1 -( 1 - 


1 ,M(N-3) 

~8~ P 4 


(17) 


(18) 

(19) 
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( 20 ) 


E 4 = T P 3 M(N ' 3) 

Although it has not been possible to obtain the general expression for the 
probability of a group containing an arbitrary number of broken elements it is 
contended that this expression is only of academic interest. It is argued that 
once the probability of secondary breaks (i. e. the probability of groups con- 
taining two. or more fractures-) becomes singif leant failure of the composite 
can be expected. More will be said of this point later on. 


Comparison with Experiments 

The ultimate test of any theory is its agreement with experimental re- 
sults. The present theory- has been compared with the data obtained in Refer- 
ences 17and 4 for glass-epoxy composites. 

The tests reported in Reference 4 were run on composites consisting of a 
single layer of 3 1/2 mil E-glass fibers embedded in two resin systems (B and 

i 

C). Tests were run to determine the strengths of the fibers for several gage 
lengths. As in Reference 4, it is assumed that the cumulative distribution 
function- can be characterized sufficiently well by the Weibull distribution. 


F(ct)'= 1 - e 


-aLj 


P 


( 21 ) 


.where L is the length of the fibers and a and j3 are parameters that character- 
ize the distribution. Using this expression the variation of'-mean fiber strength' 
with gage length is 

ct t = (aL) P ru+l/p). (22) 

.Ll 
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The fiber parameters are found to be 

- JL 

P = 9 . 40 ,* a g = 181. 5 ksi. 

, ■ From the photographs of the specimens under polarized light the inef- 
fective lengths for series B and C were found to be 0. 031" and 0. 086", res- 
pectively. 

Using this data the expected number of single broken elements were 
calculated from Equation 4, and are represented by the dashed lines in 
Figures 9 and 10. The number of fractures observed experimentally in the 
various specimens in the two test series are presented for comparison. It can 
be seen that for low stress levels there are generally more fractures than the 
theory predicts.. However, the behavior is most important at the higher stress 
levels in the area of failure loads, and here agreement is fairly good consider- 
ing the statistical spread in fiber properties and the experimental uncertainties 
involved. ' The relatively large number of fractures^at low stress levels is 
.possibly a result of damage to the fibers during fabrication of the specimens 
since glass is notoriously sensitive to handling as far as strength is concerned. 
However, since the observed fractures approach the expected number of higher 
stress levels they are of little importance. If, on the other hand, failure 
occurs without a significant accumulation of fractures the breaks at low stress 
level may be of extreme importance. Mor.e will be said of this later on. 

The quantities of E^, E^, E^ and E^ (where E. represents the expected 
number of groups of fractures having at least i broken fibers) obtained in 
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Equations 4, 13, 17 and 20 are plotted for test series B and C in Figures 11 and 
12* It can be seen that the expected number of multiple fracture groups (E^, 

E^ and E^) rise sharply in the observed range of composite failure. The fail- 
ure predictions of Reference 3 are presented for comparison. 

In order to assess the validity of the expressions for multiple fractures 
the number of groups of multiple fractures were counted on films of the tests. 
The results are plotted in Figures 13 and 14. The dashed curves in the figures 
are calculated values of E^. It can be seen that, in general, multiple breaks 
begin to appear in the stress range predicted by the theory. Furthermore, 
the composites fail without the occurrence of a’ large number of multiple- 
break groups compared to the number of isolated single fractures. 

There are several factors that could account for this phenomenon. One 
possibility is illustrated by Figure 15. This graph shows the relative behavior 

of the sum p . + p . for static and dynamic over stress factors. It will be 
2/1 of 1 

recalled that this sum represents the probability of th e 'f r a ct u r i n'g"of &t least' 

one element adjacent to a "single broken fiber/ i n t h e" ofes erve d ^ ai f uire C r a hg e 

t>jA? x j. r j'i. no6r fracture i® to be 

the dynamic curve is markedly higher than the static curve indicating that 

• .there is a definite possibility of a failure crack being caused by the dynamic 

c 

effects of fiber fracture rather than by cumulative static probability of 
failure. It should also be noted that the use of overstress -factors is just an 
approximation and that the actual stress concentrations caused by multiple 
breaks may be, and probably are, much more severe than those of single 
fractures. [As an analogy, the stress concentration factor for anellipse 
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increases with the aspect ratio) Whatever the mechanism, experimental ob- 
servations seem to indicate' that multiple breaks tend to occur shortly before 
composite failure. 

The present theory was also compared with two tests on continuous 
• (17) 

glass fibers in epoxy r.un by Friedman . The Weibull parameters for the 
strength distribution of fibers used are 


* - JL 

0 = 4. 0 , a P = 137.1 . 

'Since f3 is' an inverse measure of dispersion, these fibers had a much wider 
spread in failure stress levels than did those’ of Reference 4, This large disper- 
sion is reflected in the wide spread of the curves of E_, EL, E._ and E A pre- 

1- 3 4 

sented in Figure lk However, failure in both specimens occurred quite 
near the stress level for which the first multiple fracture is predicted. 

On the basis of the experimental evidence cited it is proposed that the 
failure stress of a continuous fibrous composite loaded intension parallel 
to the fibers can be reasonably well predicted by that load for which the first 
multiple fracture is expected to occur. That is, 


E 2 =^rV)(P 2 /i + P 3 / 1 ) = l.. 


(23) 


Analysis of Non -Cumulative Fracture Mode 
Although a large number of isolated, fractures are observed in glass - 
epoxy composites, this is not the case for other fiber-resin systems such as 
boron-epoxy and boron aluminum. For these composites failure-usually 


45 - 



(occurs catastrophically, without an accumulation of isolated fractures. How- 
(19) 

fever, Lenoe has demonstrated scattered fractures in a boron-aluminum 
-fomposite containing 5% volume fraction of fibers, which is quite small. 

‘S • 

The absence of a significant number of isolated fractures seems to indi- 
cate that the entire composite is failing at the load at which the weakest 
fiber breaks. This would mean that the matrix is actually detrimental since 
the strength of a bundle of fibers, without a matrix, would be stronger. It 
was, therefore, decided to determine the theoretical value at which a first 
fiber fracture is expected and compare the results with experimental evidence. 

Consider a population of fibers of length'!, whose strength is character- 
ized by the probability density g(a). For a sample of N‘ fibers from this group 
the distribution functioti for the strength of the weakest fiber has the following ' 
form 


(a) = Ng(cr) [l-G(a)] 


N-l 


• (24) 


Assuming a Weibull distribution for G(ct) the expected value (mode) for 
the first fiber fracture is found to be 


1 

-1 \P 


NLa (3 


(25) 


(17) 


Friedman ran two tests on boron-epoxy specimens E inches long which 
contained about 90 fibers. The Weibull parameters for the fibers were found 
to be 
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1 


f = 7. 0, a P = 368. 

The fiber stress levels, at 'failure for the test specimens were 193 and 215 ksi 
The expected value for the first fracture computed from Equation 26 is 171. 5 
ksi which is about 16% lower than the average failure stress of the two speci- 
mens. -On the other hand, the first multiple fracture is predicted (E^ = 1) to 
occur at 300 ksi while the theory of Reference 3 predicts a failure stress of 
428 ksi. In the last two cases Friedman's definition of ineffective length was 


used. 

(2 0) 

Grinius also ran tests on boron fibers in an epoxy matrix. These 
specimens were 2. 5 inches long and contained 25 fibers. The Weibull para- 
meters for the fibers were found to be 
1 

p = 11.11, a P = 433. 

Unfortunately only one undamaged specimen was tested. This specimen 
failed at 304 ksi fiber stress. The expected value of stress for the first frac- 
ture is 296 ksi while the first multiple fracture is predicted at 330 ksi and the 
failure stress predicted by the theory of Reference 3 is 340 ksi. It should be 
noted that to obtain the last two values, the definition of the ineffective length 
presented in Reference 3 was used. 

From the experimental observation of the absence of cumulative fractures 
and the good correlation between the observed failure stress levels and 
■ those predicted' for the first fiber break a good case can be made for the 
hypothesis that composites exhibiting this type of failure are only as strong 
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in tension as their weakest fiber. If this is also the case for the three dimen- 
> 

I 

s'ional composites, the matrix'must be designed so that it will prevent this type 

t • * 

of failure from occurring. This can possibly be accomplished by using a ductile 
matrix or by allowing the matrix to partially debond from a broken fiber. 

Conclusions 

A statistical model including the effects of stress concentrations for planar 
arrays of fibers in a matrix has been presented which provides a good descrip- 
tion of composite behavior up to the failure load. The stress level for which the 
first multiple fracture is expected to occur has been proposed as a predicted fail- 
ure stress. 

The model predicts that the composite itself is a "brittle" material in that 

\ -N 

its strength decreases as the length' or width of the specimen increases. This is 
illustrated in Figure 17 where the variation of with length is presented for the 
composite system used in test series B of Reference 4. This is in contrast with 
the theory in that paper which predicts a composite strength that is independent 
of length for large values of M. 

Finally, it has been shown that the failure stress level in tests of- three com- 
posites that did not exhibit cumulative damage qccurred at about the expected value 
of stress at which the first fracture was predicted. This failure stress is lower 
than that for a bundle of fibers of the same length and number indicating that the 
matrix may have a detrimental effect in composites exhibiting this type of failure. 

The present study demonstrates that the understanding of composite be- 
havior is a continually evolving process and that even in the case .of such a simple 
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loading condition as pure tension the failure mechanisms are not completely 
understood. More work must be done, for example, to explain the apparent 
difference in failure modes observed in glass and boron fibers, and to extend 
the analysis of stress concentration effects to three-dimensional composites. 
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'I'he Steady-State Viscoelastic Response 

A systematic study'of viscoelastic behavior of fiber -reinforced compo- 
sites was initiated by Hashin (Ref. 21). Such analyses are motivated by both 
theoretical and practical points of view. Fibrous composites consisting of 
linear viscoelastic phases are a natural extension of the linear elastic model 
in which phase materials are assumed .to be Hookean Solids, On the other hand, 
phase materials in composites do exhibit very strong time dependent properties 
especially in a high temperature environment. A basic understanding of the 
viscoelastic behavior of composites is important before such materials are 
utilized in practical applications. 

Hashin ('Ref. 21) related the effective viscoelastic properties of a compo- 
site which is composed of -linear viscoelastic phases to those of its constituents 
by a correspondence principle. This same model used by Hashin is used to 
study the steady state response to some simple types of sinusoidal surface 
tractions or surface displacements. It should be noted that in Ha sin's approach 
to viscoelastic theory of composites, inertia forces are- neglected so that only 
quasi-static motion is treated. The present study is based on Hashin's work 
(Ref, 21) and the same assumption is implicit so that wave propagation phenome 
na will not be considered here. . ' 

Following Reference 21, the general macroscopic viscoelastic behavior of 
a composite can be described by the following constitutive equations: 
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or 


yt) = 


r •? ic'w 

J .... (t-r) — : 
ljk'u dr 


d W T) dr 


where cr..(t) and e..(t) are, respectively, the average stress and average strain 
ij ij 

% ❖ 

in the composite at time t; G_^(t) and are the effective relaxation 

moduli and effective creep compliances, respectively. 

By the symmetry of the average stress and average strain tensors, the 
following symmetry relations hold: 




J jikt (t) ~ J ijk-t {t) “ J ij^k (t) 

In order to establish the correspondence principle between viscoelasticity 
and elasticity, it is assumed in Reference 21 

G iW (t) = G k-tij^ 


V" )= w t) 


which .are the Oasager Reciprocal Relations {Ref. 22) 

\ 

The one-sided Laplace transform of (1) and (2) gives 


V p) = p G iM <p) V (p> 
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( 6 ) 


.Jij (p) = p ^ju {p) V p) 

where p is the transform variable and the circumflex " A 11 above a function 

denotes its Laplace Transform which is assumed to exist* Because of the 

formal resemblance of* (5) and (6) to the generalized Hooke's Law in elasticity, 

A . ' . 

^ ^ijk'k^ an< * ^ijk£^ are termed transform domain (TD) effective moduli and 
compliances, respectively. 

Now for unidirectional fiber-reinforced composites, transverse isotropy 

will be considered (Ref. 21, 11}. Consequently, only five effective relaxation moduli 

(or five effective creep compliances) are independent so that (5) and (6) can be 

much simplified. For example, (5) can be written down in ter-ms of the following 

five independent TD effective relaxation moduli: 
a# 

p K^te) the plane strain TD effective bulk modulus; , 

• p G^(p) ~~ ■ the TD effective transverse shear modulus; 

P (p) -- the TD effective in-plane shear modulus; 

A * ■ 

p E* (p) -- - the TD effective Young's modulus; 

A# * ’ 

■ ^ ^ll^"*” to associated with uniaxial stress in fiber direction 

with transverse deformation prevented- by a rigid 
enclosure. 

On the other hand, the TD effective creep compliances are the reciprocals 

* 

of the corresponding TD effective relaxation moduli. For example 
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(p) 


p k,; ( P ) 


/\ * 
Pg 23 


(P) = 


P G Z3 ^ 


etc. . 

After the effective characteristic functions (relaxation moduli and creep 
compliances) are defined, Ha shin (Ref. 21) used a correspondence principle to relate 
the effective viscoelastic characteristic functions to the effective elastic moduli 
of a duplicate composite body with elastic phases. 

On the other hand, the effective elastic moduli of unidirectional fiber- 
reinforced composites have been obtained by Hashin and Ro.sen (Ref. 11). Closed-form 
expressions for four effective elastic moduli and bounds for the fifth are ob- 
tained for the "random array" geometry. In applying the corresponding prin- 

j j 

ciple, Hashin (Ref. 21) used these expressions and went through a replacement 
£c 

scheme to obtain the Raplace Transforms of the corresponding effective visco- 

elastic characteristic functions. Therefore, the entire problem of finding the 
» * ' 

quasi- static viscoelastic response reduces to Laplace Transform inversion 
which is not always easy. However, for cases where inversion is formidable, 
Abel-Tauber theorems can be used to draw important conclusions on the be- 
havior of the effective moduli and compliances without the operation of inversion. 
Refs. 21 and 11 are referred for details. 

Now, for a unidirectional fiber-reinforced composite body under special 
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boundary displacement or traction conditions (R ef. 21), (1) and (2) can be reduced to 

5{t) = f dT (7) 

J o_ 

«{t) = f* J*(t? T) — ^ * (8) ■ 

J o_ 

where a (t) is a component' of the stress (strain) tensor cr_ (t) 

whe r e a s (t) ^ J "' ( t) ^ represents one of the five effective relaxation moduli 

(effective creep compliances) of physical importance as defined previously. • 
Henceforth, (7) (or (8) ) will be used as the representative average stress-strain 
relation that defines a particular viscoelastic characteristic function. 

In the course of obtaining the steady-state response by making use specifii 
cally of the results of Ref. 21, we shall first develop certain relationships and' 
investigate their validity. Then the general results will be applied specifically 
to the unidirectional fiber-reinforced composite to obtain explicit expressions 
which characterize the response. Since we are interested in steady-state re- 
sponse to sinusoidal input, we can assume that the input'has been applied on the 
body for an indefinitely long time and that all initial disturbances have died out. 
Under this circumstance, it is' convenient to put 'the beginning, of motion at 
time -«>(Ref. 23). Hence (7) and (8) are modified to the following forms: 

5{t)= j G*(t-r) (9) - 
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and 


e (t),= f J (t-T) 


dff (T) dT 

dr. 


( 10 ) 


By changing the integration variable from T to 4 where t - T = 4 , (9) and 
(10) become, respectively, 


a(t) 


G (4) S' (t-4)d4 


(ii) 


6 (t) = j J"(4)5 (t-4) d4 

o 


( 12 ) 


where prime denotes differentiation with respect to the argument of the function. 
Using complex representation for sinusoidal oscillation, we put 

icot 


e (t) = e e 


(13) 


where is in general a complex number and to, a real number, is the angular 
frequency of the oscillation. 

Substituting (13) into' (11) we obtain 


- . iCOt 

cr (t) = icoe e G (co) 


(14) 


•where 


G*(W) J G*(4) e d4 

o 


(15) 
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is the one-sided Fourier Transform of G (t) if it exists. 

!' Let ff(t) = <J e where cr is in general a complex number, then by (14), 


g(t) 
e (t) 


i a) G*(o) ) M - " 


( 16 ) 


•where M"' is called the effective complex modulus associated with G (t)'. 
Equation (16) can also be written as 

- , , i „ i (Wt + 6 ) 

a (t) = | M | e e 


(17) 


where 6 - tan 


-1 


Im [M ] 
Re [M”] 


which reveals the fact that 5 (t) is not in phase with e (t) 


— ico t" 

Similarly, if CT (t) = cr^e then (12) becomes 

, . . ftotC*, . 

e (t) = lcocr e J (w) 
o 


(18) 


where 


r>) 


r (4)e -' it0 W . 


(19) 


is the one-sided Fourier Transform of J" 1 ' (t) if it exists. 

Let e(t) = e e i£l)t then by (18), 
o 

£ (t) = = i 60 J*(w) s (20) 

a (t) c 0 ’ 
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or 


: 41 .<*» 

which again shows that £ (t) lags behind a (t) with a phase lag 6 . 

From both (16) and (20), it follows that 

? ' /—it. 

- W G (W) J (CO) = 1 

if both G" (w) and J*(W) exist. 

The physical significance of Im[ M ] andR e [M*'] is now discussed. 

First of all, it can be shown that the rate of work done by surface traction at time t 
on a body can be expressed in terms of an integral over the entire body as follows: 


j T. u. d S = J cr.. e.. 

ij i] 


d V 


where S denotes the boundary surface and V, the volume of the body. The above 
equality is valid only under the assumption of quasi- static motion. 

Furthermore, under special boundary conditions on the boundary surface 
S (displacement boundary condition u.(S) = e (t) x or traction boundary condition 
Ti<S) = cr.^t) n , (R ef . 21)), it is easy to show that:- 


CT.. (x, t) s.. (x, t) d V = CT..(t) e.. (t) V 

y - ij - iJ y 


Therefore, the total work done from time t^ to time 1 2 . is: 


r. u. dS dt = V \ cr..(t) e .. (t) dt 

ii-J 1 iJ y 


*1 S 
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If x - is the period of motion, then 


W = V , 


2 - 

p 


^ e..(t) dt 

y 


(22) 


denotes the total work done on the system in a cycle. According to the First 
Law of Thermodynamics, for isothermal steady state deformation, W is the 
total energy dissipated and transferred to the surroundings in the form of heat 
in a cycle. Under the special boundary conditions by which (7) and (8) are ob- 
tained, W in (22) can be reduced to the following simple form: 


W = VF 


cr (t) e (t) dt 


(23) 


■ where the factor F is either 1 or 2 depending on the boundary condition. 

Therefore, if e (t) = Re [ e^e 1 ^ ] and ff(t) = R^ [M’' e^e 110 * ] , then for 
a cycle with angular frequency of motion to, after some manipulation, (23) 
becomes 27T 

r 00 . 2 

W = VF i 5 s dt = I t .Im [ M ] VF (e is assumed 

J • ' o o 

T 

h 

to be real without loss of generality ) where T is any time during the motion. 

\ 

This gives the amount of energy dissipated in a cycle. 

Furthermore, it can be shown that a fraction of the amount of energy 
S o (& e [ M"‘]j VF is twice in. every cycle alternately stored and expended in 
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the system. Therefore, Im [M ^ ] is related to the energy dissipated while 
Re related to the energy stored. 

From the above results, it is clear that in the steady state, the 
responding average stress (or average strain) will vary sinusoidally with the 
same angular frequency 0) as that of the input average strain (or average stress' 

with a phase difference 6 which is a function of 0 ) . Moreover, the amplitude 

u • 

ratio is also a function of CO only, 
s 

o 

According to (16) and (20), the knowledge of hinges on the knowledge ’ 

‘ r\-. 

of G (to) or J^(co) . However, G n (C0) and J (CO) are formally related 

/ V 

to the Laplace Transforms G P (p) and J ''(p) by. the following' relations: 

, j 

G* (CO) = G*(i(d) (24) 

and 

a! 

J v (co) = j"(i CO) ' (25) 

/'.I. ^ .1. ' 

Therefore, if expressions for G (p) and J (p) are known, G (CO) and J (to) 

can be obtained immediately by the replacement of variable p' in G' (p) and 

A,. ... 

J (p) by iW , However, there~s till exists the problem of convergence -- 

given a function of time t, the existence of its Laplace Transform does not im- 

( . 

ply the existence of its Fourier Transform. Thus, given-.a Laplace Transform 
of a viscoelastic characteristic function, it is necessary to examine the location 
of its singularities in a complex plane of complex variable p. According to the 
definitions of the one-sided Laplace and Fourier Transforms, it is observed 
that if the singularities. of the Laplace Transform are all located in a region to 



the left of imaginary axis of the complex plane, j en the corresponding Fourier 

\ 

Transform exists. \ 

•(' 

Hashin (Ref. 21} has obtained closed-form ejoressions for some effective 
relaxation moduli and creep compliances for the 'random array" model. Then, 
by making use of (24) and (25), the corresponding effective complex moduli can 
be obtained in a straight-forward manner. For example, for elastic fibers 
embedded in a viscoelastic matrix, the effective relaxation moduli (t) 
characterizing the plane strain dilatation has the following form in the Laplace 
Transform domain: 


p k 23 (p) = [k (p) + -r_( P ) ] + v 


3 T 1 K f + i G f - k m (p) - ir m (p) T k m (p) + | r m (p) 


-l 


where k^(p) , r^fp) are matrix T0 moduli; and are fiber elastic bulk 

and shear moduli; v, and v are the fiber and matrix volume fractions, respec- 
f m 


tively. All these quantities are considered known if the composite body is 
given. Then according to (24), formally we have: 

1 , ■ 


icoK (u>) = [k (io>) + — T (i«) ] + v 
23 m 3 m 


f t K f + I G f- k m (i “> - I r m (i “> k m (i “> + ! ^ 


and the complex moduli associated with K 23 (t) is, according to (16) , 

.V - — 

• Ml' = iWK (w) 

23 23 ' 


Furthermore, it follows from Ref. 21 that for rigid fibers, 


K 23 (W)= K m (W) + (“) 


f 


• ) +[K m (W) + | G m (W> ] (1 - v f ) 
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where K (to) and G (to) are the one-sided Fourier Transforms of the matrix 
m m 

bulk and shear relaxation moduli, respectively. On the other hand, if fibers 

are rigid and the matrix is elastic in dilatation and Maxwellian in shear, - 

Hashing ef.21)has obtained simple expression for the effective creep compliance 
% 

^2 3^) w kich the one-sided Fourier Transform does not exist. Similarly, 
other steady-state responses such as in-plane shear, etc. can be obtained in a 
straight-forward manner. 

It is emphasized her.e that the method developed here is also valid for 
other models than the composite-cylinder-assemblage model used in Ref. 21; the 
latter is only a case in which bounds of some of the characteristic functions in 
transform domain coincide to yield closed-form expressions. If other ex- 
pressions can be obtained by other models, the corresponding effective com- 
plex moduli can be obtained by the replacement scheme defined by (24) and (25). 
On the other hand, vibration experiments can be performed on composites to 
determine experimentally their respective effective complex moduli from ■ 
which information on their corresponding effective viscoelastic characteristic 
functions can be obtained through Fourier transformations. 
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Evaluation of Transverse Effectiveness Factors for Use in Elastic Analysis 


of Three Dimensional Filamentary Composites 

The method of analysis developed in Reference 1 for the elastic constants 
of composites having filamentary reinforcements at various angles to the three 
principal orthogonal axes of the material employed factors P 


etc. to define the transverse effectiveness of the filaments for resisting stretching, 
shearing, etc. For simplicity the assumption was made that the transverse 
effectivenesses were the same for all strains due to stretchings regardless of 
whether the strains were' -induced by Poisson's ratio- effects, or were simply 
those in the direction of an applied extension. (Different values of transverse 
effectivenesses in shear from those in stretching were, however, allowed for.) 
This approximation led to generally satisfactory results for stiffnesses, with 
the possible exception of E^» that -is the stiffness along the filaments, for uni- 
directional reinforcement - a somewhat disturbing exception in view of the 
simplicity of calculation of this stiffness by the rule of mixtures. Even more 
disturbing was the fact that the values of Poisson's ratios calculated with this 
simplifying assumption of equal transverse effectiveness factors were not con- 
sistent with those found by other methods of calculation. c 
( Guidelines for the extension of the analysis of Reference 1 to provide for 
. various values of for direct and Poisson strains were included with the 

analysis, but within the time available no evaluation could be made to determine 
whether or not the extension would be justified. A brief series of calculations 
has now been made to evaluate the differences among the various transverse 


effectivenesses for the extended analysis. The results are shown in Figure 18, 
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and the equations relating thereto are reproduced as Tables 1 to 3. 


In Figure 19,curves of 




■v 


3 LT 


, and 


are plotted as 


derived from the upper limit values of the elastic constants of Reference ]1 for 
a typical glass/ epoxy combination employing the assumption that 

e p 0 ~jP. 0 ~ / • As can be seen, differences among ^Oj,’ 

’ an< ^ fi o are ^ oun< ^> °i increasing relative magnitude with in- 

creasing volume fraction of binder (the abscissa on the Figure). Also 

differs from the value unity; this is the variation which is to be anticipated • 

from the previously noted inconsistent calculation of E for uni- directional 

J-J 

reinforcement' with the simplified transverse effectiveness assumption. The 
differences in the various ^ - values shown are, of course, just those consistent 
with the elastic constants as found from Reference 11. The .use of the p jf of 
Tables 2 and 3 thus reproduce properly all the elastic constants 
of a unidirectional reinforcement configuration, 'and hence provide a 
self-consistent basis from which the effects of multi-angular reinforcement, 
maybe determined via extensions to equations like those of Reference 1. 

These extensions have been made and the general equations for the compliances 

5 •" 

incorporating the various '^S<S are presente<hdn-Tabler4.cO.,,.. *"o. 

As presently writtqnj^the ,^qu%tip.ns ;j b£.T^bJftS'.2^ 3, aiidh-4 conceal- as ome 
internal inter- ; r.e^jppshjps n ^mqpg ; the various; expressions, aFiOii.example, the 
trigonometric expressions of Table 4 are for convenience written in terms of 
the three angles ^ , fjj? , and jfL which the reinf.orcing filaments make 
with the composite orthogonal axes of symmetry (1, 2, and 3). Only two 
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"kngles are needed, however, to define the filamentary directions; 
.and jf]L are related by the well known formula 




a) 


Similarly, for convenience six transverse effectivenesses have been 


defined ( 


0 > 


- see tables 1, 2, and 3), 


where only five are needed for consistency with the five elastic constant: 

used in' evaluating the transverse effectivenesses. Accordingly the equations 
of Table 3 may be combined, as for example to yield the relation 



where = the Poisson's ratio of the filamentary material. 

Thus simplifications of the equations as given in Table 4 are undoubtedly 
possible by the employment of the trigonometric expressions, and ^ values 
which lead to the least complex algebra. 
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DESIGN CONCEPTS 


- Study of Composite Structures 

Efficient application of composite materials to aerospace structures 
requires proper selection of reinforcement pattern and material as well as 
; overall structural arrangement. To a large degree guidelines for optimum 
‘ design of such composite structures have been lacking, nor have they been ob- 
.vious a priori to the designer. 

Some clues about possible directions toward efficient configuration have 
recently become available. Reference 24 explored effects of material and rein- 
forcement pattern on the structural efficiency of boost vehicle shells. Reference 
2 investigated similar effects for a number of aircraft structural elements. 

While these studies revealed a number of specific factors of importance for 
efficiency of application of composites, three general conclusions also evolved 
which served as a basis for the studies to be reported here. These conclusions 
were: 

f 

1. For a wide range of shell type applications the isotropic (0°, ±60°) 
reinforcement configuration is most efficient. 

2. The high .{multi- directional) stiffness-to- density ratio of beryllium 
makes it more efficient than most near-future composites for many 

- S' 

shell and plate type structures, - particularly if the load intensities 
encountered are low, or if stiffness requirements are important. 

3. Filamentary composites appear most attractive when used as 
unidirectionally reinforced elements to carry unidirectional loads. 
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The further question that seems to arise from these conclusions is 
whether some combination of the biaxial properties ,of beryllium (or isotropic 
cpmposite) with the uniaxial properties of one-directional filamentary rein- 
forcement may possibly be the best configuration of all. 

Circular Tube -Columns 

To explore this possibility on an orderly basis, a simple round-tube 
column was selected as a first model for study. This model was chosen for 
simplicity; for an isotropic material optimization of the thin-walled tube in 
compression has perhaps achieved the status of being a classic example of 
balancing proportions between local buckling and column bending to achieve 
minimum weight. ' Thus its use to explore effects of combinations of unidirectional 
and multidirectional properties is uncomplicated by complexities of the optimiza- 
tion procedure itself. 

Accordingly beginning with the model of Reference 25, modifications were 
assumed of increasing complexity, and the effects on efficiency evaluated. As. 
a first step, unidirectional'ly reinforced composite stiffeners were assumed added 
to the tube at three equally spaced points around the circumference. These 
stiffeners were assumed to contribute nothing to the local stability, -or instability,- 
of the tube walls, but were' assumed to add effective areas for column bending 
so that the effective radius of gyration of the tube is given by the expression 
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, f 'radius of gyration of stiffened tube 
outside radius of tube 
thickness of tube 

Young's modulus of unidirectional (boron) composite 

Young’s modulus of tube material 

volume fraction of composite stiffeners 

volume fraction of tube material (V D + V = 1) 

B o 


With equation (1) incorporated into the procedure of Reference 25, the 
efficiencies of the reinforced tubes were calculated from the formulas that 




Gyr» Eff 


R 

t 

E„ 
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axial load 

v* 

. column length 
overall weight 

density of stiffening material 
density of tube material 

and the subscripts Sec and Tan refer to the secant and tangent moduli of the 
tube material, respectively. 

As indicated in the formulas, the reduced moduli used for plastic buckling 
were the tangent modulus for column bending, and the secant modulus for local 
buckling. The use of the tangent modulus for column bending is well founded. 
The use of the secant modulus for local buckling is used here as slightly more 
optimistic than the root mean square of the tangent and secant moduli as pro- 
posed in Reference 26. Thus the efficiency curves resulting are perhaps slightly 
too high in the plastic region, representing a kind of upper bound. (Most 
affected are the curves for beryllium which should be accordingly somewhat 
discounted at the upper end. ) In any event the 'use of this possibly optimistic 
reduced buckling modulus changes none of the conclusions drawn from the 
results, 

Results of the calculations are given in Figure's 19 to 24. In Figure 19 are 
given basic results for unstiffened aluminum-alloy and beryllium circular-tube 
columns. Here.the beryllium tubing is found to vary from 25% of the weight of 
the aluminum- alloy tubing in the elastic range 



whpe 
■ P 
' L 


W 

s 
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to 80% of the aluminum- alloy weight at high stresses 


kM 

The addition of the boron/epoxy three-point reinforcement permits the 




aluminum- alloy tubing to be made substantially lighter both at high and low- 

loadings (Fig, 20).- In the elastic stress range there is apparently an optimum 

/ % \ 

reinforcement ratio j •**=***■ s beyond which higher percentages of 

* . 

reinforcement do not further increase the efficiency. At the high stresses, 
however, the very great strength and stiffness of the unidirectional reinforced 
composite provide increases in efficiency up to the maximum reinforcement 

/ % . i ■ 

ratio considered ff *«==» -g> ^ r 

t • 

The beryllium round tubing is not improved by the boron composite rein- 
forcement at low stresses; rather the efficiency is decreased as .the reinforce- 
ment ratio increases. This trend is just barely reversed, as might be expected, 
at high loadings (see Fig. 21). 

Comparisons between Figures 20 and 21 reveal that at low loadings the 
unreinforced beryllium is always the lightest of the combinations considered, - 
being approximately one-third the weight of the best boron/epoxy- aluminum. 

At the high stresses on the other hand the reinforced aluminum is better than 
any beryllium tube, being about one-half the weight of the beryllium construc- 
tion at the maximum loading considered. Thus it appears that the low propor- 
tional limit strain of the. beryllium together with the flat top to the beryllium 
stress -strain curve beyond the proportional limit prevents it from acting 
effectively in conjunction with high-strength uni-directional reinforcement. 

In this range a better combination of materials is an isotropic boron/epoxy 
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configuration for the tube walls. This configuration should still be elastic at 
tf 1 ® values of P / L 2 covered in Figures 19 _ 21, and at the maximum values of 

J V* 

S/juZ considered the isotropic boron/epoxy tube is calculated from equations 
(1) and (2) to weigh very nearly one-half as much as the best reinforced 
aluminum tube. 
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Scalloped. Tube Columns 


One of the pitfalls of efficiency studies like the foregoing is the possibility 
that a poor geometrical shape has been chosen, and as a result misleading con- 
clusions are derived about the effects of various materials of construction. 
Ideally an optimum shape should be used (if need be the optimum for each 
material combination) so that shape effects maybe divorced from material 
effects in the evaluations of the calculations . 

In order to assess the shape effect, the efficiencies of scalloped tubes 
were also investigated. First, a three -lobed scalloped round shape was used 
(Figure 22). The objectives of the scalloped shape are enhanced local buckling 
resistance due to decreased radius of curvature of the tube wall together with 
increased column bending resistance because the scalloped .tube approximates 
an equilateral triangle, and' the equilateral triangle has a 21% greater cross- 

sectional moment of inertia for the same area than a circle. Because the 

> 

scalloped round tube does not substantially enhance the moment of inertia, . 
however, it was abandoned in favor of a scalloped triangular tube. The charac- 
teristics of these two shapes are summarized in Figure 22. 

In Figure 22 are plotted the percentage increases from a simple thin-walled 
■ ' c . 

round tube section attained by scalloped-rOund and scalloped -triangle sections 
having the same cross-sectional area and radius of curvature-to-thickness 
ratios as the reference round tube. As indicated by the curves of the Figure, 
the "scalloped triangle" has substantially greater potential for enhancing the 
column strength (at a constant so that in first approximation the local 
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buckling strength is constant). than the "scalloped circle." Furthermore, the 

apexes o£ the triangular shape provide a greater radius of gyration for added 

r- 

stiffening as in the form of concentrated boron/ epoxy elements than the crests 
of the scallops of the "scalloped circle. " 


• ‘While clearly the extreme limit of 21% increase in section properties 
represented by the "infinitely, scalloped" triangle ® @j! °£ Figure 4 surely 

cannot be attained, the development of a reasonable fraction of that amount may 
be anticipated for "reasonable" proportions. Here "reasonable" is taken to 
mean that the cusps between scallops are of sufficient- depth to establish stable 
corners not prone' to local buckling at the stresses encountered by the tubes. 

The assumption was made that the proportions so labelled on Figure 22 are 
reasonable, and likely performance gains through their use for reinforced com- 
posite columns was calculated with the results given in Figures 23 and 24.- For 
these calculations, equation ( 2) was replaced by 



where 


moment of inertia' of scalloped triangle section 
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area of triangle section (equal to area of circular 

section of equivalent } 

inside radius of scallops 

equivalent thickness -radius ratio of tube 

. /,... 

In Figures 23 and 24 comparisons- of the efficiencies of the reinforced 
scalloped triangle tubes with envelope curves representing the lightest reinforced 
round tubes (the dotted curves on Figures 23 and 24, derived from Figures 20 and 21) 
show the expected gains for the triangular shape. With the triangular shape 

higher reinforcement ratios for the aluminum alloy tubing are effective (Fig. 23) 

% 

and a reinforcement ratio «===» = 4 produces a tube column competitive with 
beryllium even in the elastic range (see Fig. 21). Further, the boron/epoxy re- 
inforcement of beryllium in the plastic range is more effective with the scalloped- 
triangle shape. The general trends established in the studies of the reinforced- 
round tube columns, however, are not changed by the mo r coefficient triangle 
shape. Thus the scalloped-triangle beryllium tube is the most efficient of all at 
low loadings and the reinforced aluminum -alloy becomes more efficient when the 
stresses would cause yielding of the beryllium. 

The gains in efficiency shown on Figures 23 and 24 for the scalloped-triangle 
shape appear great enough so that experiments s to determine their validity appear 
desirable. That is, experimental definition is needed of the stability of a 
sc all oped -wall tube to establish the degree of scalloping beyond which corners 
of cusps will no longer remain fixed as straight lines along the tube to provide 
edges as restraints against local buckling distortions.. 
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BTfrte Efficiencies 

j 

; The efficiencies of composite materials for plate applications, already 
investigated for a variety of material combinations in Reference 25 , have 
exhibited similar characteristics to those already discussed above for tube- 
columns. Indeed combinations of beryllium and 0° boron/ epoxy plates made 
up as sandwiches with the beryllium as the faces and the 0° boron/ epoxy as 
the core demonstrate rather clearly the problems encountered in the develop- 
ment of combination metal composite structures, as is shown in Figure 25. 

In Figure are plotted the results of calculations (made as in Reference 

25) of the efficiencies ^of sandwich plates which have 0° reinforced boron/ epoxy 
cores and faces of beryllium, isotropic boron/epoxy, and (at the very top end of 
the curve) no faces at all. The solid curve represents the beryllium -faced 
sandwiches; up to the discontinuity (at J ) tiie plates 

are totally elastic (to achieve this elasticity at such high values of a com- 

, ~~b " 

pres'sive pre -stress is assumed in the boron/ epoxy core just sufficient to stress 
the faces to the proportional limit in tension at zero external load). The rapid 
increase in weight of sandwich above the discontinuity arises from the continuing 
decrease in stiffness of the beryllium faces as they are stressed further and 
further beyond their proportional limit. The dotted curve represents the higher 
efficiency that would be attained with beryllium with a higher elastic limit (or 
pre -stress). Inasmuch as neither of these improvements in the beryllium 
response appear possible of attainment in practice, the solid curve is the 
realistic one for the sandwich, and, as indicated, it can be surpassed in 
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efficiency — as by sandwiches utilizing isotropic boron/epoxy faces at the higher 
load intensities. Indeed as the loads become high enough, only the 0° boron/ 
epoxy retains enough stiffness to be effective, and no faces at all become most 
efficient. 

To a degree, the curves of Figure 25 may be considered to depict in 
general the y^rious. characteristics which must be balanced in composite plate 
construction for maximum efficiency, as follows: 

(1) 0° reinforced material should be used as the strength element, 
if possible pre -stressed to provide the maximum possible 
elastic range to the material used to provide transverse 
stiffness or plate buckling resistance. 

(2) The material used to provide the plate buckling resistance 
(the facp material in a sandwich) should be selected on' the 
basis of the stresses to be carried, — at’ low stresses- beryl- 
lium is most attractive, at higher stresses composite con- 
figurations approaching closer and closer to the‘'0° config- 
uration should be used. 


(3) While the example shown (Figure Zfyis for a solid sandwich, 

M? fiL fjk | 

so that the average stress is simply s S , the 

t bit# 

same rules as (1) and (2) apply for hollowed out constructions 


which raise the average stress to ^x where t is the aver- 

1 

age. thickness and t t . Thus, in general, any hollowing 


out (lightweight core sandwich construction, etc.) which 
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does not introduce new instabilities or failure modes is 
like an increase in the structural index value on Figure 25, 

2r 

- i.e. , it tends to make optimum constructions utilizing 
higher percentages of 0° reinforcement and low.er 
percentages of beryllium. 


•** ' efiiciehcie y OI * — B/ c. .id boryl.J#*;i-okii Z-p 0 i si JV.oti ucvio 


These comp. risons are presented in Fi<r. 71 
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Panel Efficiencies 


The plate efficiencies just discussed were considered on the basis of 

initial buckling stresses only. Particularly at low stresses, as is well known, 

plates of most materials can continue to carry compression load without failing 

jit stresses substantially above the initial buckling stress. 

< 

For plate assemblies incorporating 0 reinforced composites as stiffeners, 

initial buckling of the plates may be expected to cause the major part of any 

subsequent load increase to be borne by these stiffeners. This load transfer 
o 

into the 0 reinforcement is precisely the mechanism noted desirable above for 
■ increased structural efficiency. Accordingly, for example, compression 
panels incorporating 0° reinforced stiffeners should be expected to achieve 
high structural efficiencies. 

High structural efficiencies were indeed found for 0° boron/epoxy Z-section 
stiffeners on ±45° reinforced boron/ epoxy skin in Reference 2. Part of this 
high efficiency arose from the high shear stiffness of the skin material 
(allowing the use of thin skins and hence allowing most of the material to be 
incorporated in the stiffeners.) 

Reviews of the properties of ±45° boron/epoxy and comparisons with the 
properties of beryllium suggest that replacement of the ±45° boron/epoxy skin 
on the panels of Reference 2 with beryllium, should result-in even higher effi- 
ciencies than those calculated therein. For example, pertinent properties of 
these two materials are compared in the following table: 
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±45°B/Epoxy 

’ Be 

_ . ' Mg 

Density — 

2. 19 

1.86 

' 

(pci) 

(0. 079) 

(0.067) 

MN 

> 


Yield Stress — =- 

138 

400 ■ 

(ksi) 

(20) 

(58) 

GN 

Shear Modulus — =- 

' 72.5 

145 

(ksi) 

(10,500) 

(21,000) 

Y oung 1 s Mo dulus GN 

25.5 

304 

2 

m 

(3700)-- 

(44, 000) 

(ksi) 



Yield Strain 

0.0054 

0.0013 


Thus, except for its low yield strain, beryllium appears vastly superior to the 
±45°B/Epoxy. The low yield strain, however, raises uncertainties, and a 
detailed evaluation is required to determine whether the beryllium/ composite 
panels are as superior as would at first appear. 

Such an evaluation was made as a part of the study of Reference 27. The 
results are reproduced here and extended to make possible direct comparisons 
of the efficiencies of the ±45° B/Epoxy and beryllium- skin Z-panel constructions. 
These comparisons are presented in Figure 26. 

s 

Figure 26 plots curves of weights of box-beam compression covers (plus 

x , ^ 

supporting ribs) for optimized Z-stiffened panels. Optimization includes the 
selection of rib spacing and stiffener size and spacing for minimum weight to 
carry a design bending moment M with a box beam of width b and depth d.- 
Further the skin thickness, is required to be adequate, in terms of the box.depth. 
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, to provide a specified torsional stiffness; - accordingly the beryllium and 

J ±45° B/Epoxy are compared at "equivalent" values of skin thickness t 

• S Eq 

v such that 



The curves of Figure 26 show substantial weight savings for the beryllium 
skin on 0° B/Epoxy stiffeners. These weight savings are depicted by the cross- 
hatching between the curves for the ±45° B/Epoxy- skin panels (the dashed 
curves) and the beryllium- skin panels (the solid curves) at the equivalent tor- 
sional stiffness measures (equal values of ). 

.... 

Several factors accumulate to produce the weight savings shown in Figure 
26 for the beryllium- skin composite construction: (1) the favorable shear stiffness 
of the beryllium permits an even higher percent of stringer material-with its 
high-strength, high column-bending stiffness characteristics, - than for the. 
±45° Boron/Epoxy skin; (2) the beryllium is of lower 'density than the boron/ 
epoxy; (3) the high strength of the 0° reinforcement (as noted in the opening dis- 
cussion for. this section) permits it to carry the load long after initial buckling 
of the beryllium, skin. 
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Studies of Composite Materials 


Three -Dimensionally Isotropic Materials 

f ■ The form of the matrix of elastic constants relating stress to strain in a 

generalized Hooke’s law has been studied for various conditions of structural 

symmetry. For two dimensions, it is known that a material is isotropic in a 

plane for which there exists a normal axis of at least six-fold symmetry. Thus, 

as a practical example, a fibrous composite plate is effectively isotropic in 

its plane when one n of the fibers are oriented every — radians for n = 3. 

n 

(By effectively isotropic,- it is meant that the average stress, average strain 
relations are isotropic. ). Similar conditions to obtain three dimensional iso- 
tropy appear to be unavailable in the literature. • Several configurations having 
multiple symmetry conditions which might yield at least an approximation to 
an elastically isotropic material have been studied. The first material treated 
has. two three -fold axes of symmetry separated by the angle for which the cosine 
is ’-1/3. Such a material can be obtained with a fibrous composite having one 
fourth of the fibers oriented in each of the four directions defined by the altitudes 
of a regular tetrahedron. 

The stress -strain relations for a general elastic body may be written as: 

r.. = c.. e . i. 0 = 1 , 2 .. .6 (1 ) 

where 7 • are the stress components 
1 

S are the strain components 

3 

c.. are the 21 independent elastic constants (c.. = c..) 

13 Ji ij 
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We treat a material having two axes of three -fold symmetry at an angle 0 
.such that • 


cos 9 = 

As a consequence of this, there must exist two additional axes of three -fold 
symmetry; the four axes being oriented along the directions of the four alti- 
tudes of a regular tetrahedron. (This four-legged array is called a caltrop.) 

i rr 

We select the vertices of the tetrahedron at the points: (0, 0, — )> 



This is a tetrahedron with centroid at the origin. This body has elastic symme- 
try with respect to a rotation of — ^ about the z axis. Thus for the following 
transformation of coordinates, the stress-strain relations remain unchanged: 
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This array has seven independent constants. 

The stress- strain relations will also be invarient for a rotation about 

one of ths other axes of three-fold rotational symmetry. In particular, we 
i 

consider a 120° (counter-clockwise) rotation about the axis passing through 
the point (jfi , -~ . - followed by a 120° (clockwise) rotation about 

' ' / / i 1 

the axis passing through the point/ yrg— , 0, - 

for this transformation are 


jl 

4*3 


The direction cosines 



X 

y 

Z 

X 1 

1 

3 

0 

2/2 

3 

y» j 

0 

-l 

0 

z‘ 

2/2 
z ' 

o j 

-L 


It can be shown that the elastic constant matrix must therefore be of the 


form 


c.. - 
3J 


11 


12 

V C 12- C 33 ° 

/2( c u-c 3 3) 

■; 0 

ll 

C U +C 12- C 33 ° 

v ' 2(c ir c 33 ) 

: 0 


o 

CO 

co 

o 

0 

0 


3 

T c ir 

=£-c o’ * 

2 33 

-/2(Cu 



'3 C 12 




2 c n" 2 " 

°33 °- 


2 (C 11'°12 ) 


This material has only three independent elastic constants and is a 
material with the symmetry of the cubic crystal. Indeed, a transformation 
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■ of the principal axes of a cubic crystal given by: 



x" | 

y" 

z" 

x' 

A 

1 

72 

0 

y' 

1 

76 

1 

7T 

-if 

z' 

1 

73 

i 

7T 

1 

73 


followed by a transformation of axes given .by - 



x' 

y' 

z' 

X 

0 

. -l 

0 

y 

1 

0 

0 

Z 

0 

0 

1 


brings the principal axes of a cubic crystal in coincidence with the x y z axes 
which are the reference axes for eqs. (5). This transforms the array of elas- 
tic moduli of a cubic crystal to the following form; 

0 -jVtC-D-ZE) 0 

' 0 “ITS < C ' D - 2E > 0 

' / 

0 0 o 

\ 

j-(C-D+E) .0 ~g.(C-D-EI 

-i(C-D+E) 0 








|i(C+D+2E) -g- (C+5D-2E) y(C+2D-2E) 


c.. - 

aj 


-(C+D+2E) — (C+2D-2E) 


— (C+2D+4E) 


-f(C-D+4E' 

o 

(4 ) 
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where C, D and E are the three independent elastic moduli for the cubic 
system. /* 

The array (4) has the same form as the array (3). - In fact, if we let 


c = 4< *n " 3c 33 


D = -c^ + 6c + 6c 0 


and 


12 33 

„ (C U + C 12 ) . 

E ~ ' 2 +C 33 


then (3) and (4) are identical. This shows that the material with the elastic 
moduli of the form of (3) is elastically cubic. 

Although this material is not isotropic, the properties for a particular 
fibrous composite prove interesting, as will be shown in the following sec- 
tion. 

Another configuration of interest is that of Reference 1 which is construe 
ted by taking three pairs of reinforcing filaments, oriented with reference to 
an orthogonal Cartesian coordinate system, xyz, as follows: one pair in 
the xy plane making angles of +0 with the x axis; one pair in the yz plane 
making angles of i-9 with the y axis; and one pair in the zx plane making 
angles of +& with the z axis. It was shown in Reference 1 that when 0 is ap- 
proximately 30°, the shear modulus, Poisson's ratio, and Young's modulus 


satisfy the relation: 



It can be shown that when 


0 = tan (2 sin> 18 ) & 31° 43' 

/ 

the twelve rays ( 6 lines) from the origin pass through the vertices of a 
regular icosahedron with centroid at the origin. These six lines are there 
fore axes of five-fold symmetry. The possibility that a body having six 
axes of five-fold symmetry.as isotropic is suggested. This result is per- 
haps of academic interest .only, as it does not appear possible to construct 
a continuous space lattice. having five-fold axes of symmetry. 
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Evaluation of Moduli of Caltrop Reinforced Materials 
The equations derived for the elastic constants of three -dimensionally 
reinforced composites were used to compute properties of the caltropic re- 
inforcement configuration. Calculations were made using both the equations 

9 . 

of Reference 1 (for which ,j3 =(3=3 , etc. } and the equations of Table 

°T ®LT • °G 

4 herein. The following material constants' were employed.' 


(Young's modulus of filaments) 


10 


V i 


(Young's modulus of binder) 

V 

/• 

(Poisson's ratio of filaments). 
(Poisson's ratio of binder) 


2T xE f 

0.2 

0. 35 


These values correspond to those for E-Glass in epoxy. Values of 3 were 

found for these material constants by derivation of the upper bound Hashin- 

/ 

0 

Rosen (Ref. 2) elastic constants of unidirectional reinforcement at values of 

volume fraction filament of 0.2, 0.4, 0.6, and 0.8. 

Results of the calculations are given in Figure 27 for the filament orien 

tation shown. Both methods of calculation yield different compliances in the 

2-and 3-directions (values of A. and A.) from that in the 1-direction. Differ 

4.-6 

_ ences between the twq methods .are not .substantial-, .as was, anticipated, .nor 
are the compliances, found.s.uhstantially. at variance in^the -thr.e.e -directions., 

F;Q £ engmeer‘ing - .p,urposes, the compliances in the three directions are pro- 
bably so nearly equal that differences among them can be neglected. 
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Three Phase Composite Compression Members 

In the foregoing studies of the efficiency of combination metal /composite 
structures the efficacy of 0° reinforced composite material to provide strength 
to the structure while, some other material or configuration provided adequate 
transverse stiffness or continuity was evaluated. On the presumption that such 
0° reinforcement would truly be an extremely high compressive strength material, 
these evaluations were indeed favorable. Because of the relative weakness and 
lack of stiffness of available resin binder systems, however, it is not obvious 
that the 0° configuration necessarily has the maximum compressive strength. 
Rather, as this' section will demonstrate, some transverse reinforcement may 
be desirable. 

Reference 4 described the mechanics of stabilization of uni -directional 
filaments in compression. As noted there, three failure modes may be encoun- 
tered: 

\ 

(1) a mode-in which the shear stiffness of the binder material governs, 
described by the equation 



where 


% 


f 


ultimate compressive stress in filament direction 
shear modulus of binder 
volume fraction filament 


(2) a mode in which the strength of the binder is critical, described by 
the equation 
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where 



yield stress for binder 
Young's modulus for filament 


and (3) a mode in which the stiffnesses of both filaments and binder are 
operative, described by the equation 



The lowest value of gu found from the foregoing equations is, of course, the 
€> " 

failure stress. 

In all three modes of failure the binder properties are important. Accord- 
ingly, increases in binder properties may be expected to be reflected in increases 
in compressive strength of uni -directionally reinforced composites. One method 
for increasing effective binder properties would appear to be to take some of the 
uni -directional, load-carrying filaments and distribute them uniformly through- 
out the binder. Thus, while for every failure model the value of would 

tend to be reduced by the reduced volume fraction of load -carrying filament 
(V£ in above equations), the net effect should be expected to be an increase in 
££ due to the improvement in binder stiffnesses and strengths,. 

In order to assess the potential magnitude of compressive strength 
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increase available by this approach, the assumption was made that the binder 

properties attained by the utilization of some fraction v of the- filamentary 

*b 

material {the original volume fraction filaments v- was assumed to be 0.5 
* . *o 

throughout) for binder reinforcement would be given by simple, rul e -of -mixture - 
like, expressions such as 



where 

^ is the Poisson's ratio found by the rule of mixtures 
for filaments and binder 
and . 

with the primes denoting the properties after. binder reinforcement. 

Using the foregoing equations, three possible composite combinations were 

evaluated: (1) boron binder reinforcement for uni -directionally reinforced boron/ 

epoxy; (2) glass binder reinforcement for uni-directionally reinforced glass/ 

epoxy; and (3) glass binder reinforcement for uni-directionally reinforced boron/ 

epoxy. In the first two cases the total amount of reinforcing material in the 

composite was held constant, so v„ + v = 0.5. In the last case v. was held 

b * - 

constant and the total amount of reinforcement therefore increased as v was 

D 

increased. 
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The results of the calculations for these three combinations are plotted 

in Figures 28, 29, and 30. In Figure 28, curves for the boron-boron com- 

■S 

posite are presented for the following nominal constituent properties: 



^5 = 0.2 

The boron- glass combination of Figure 30 uses the same nominal properties 
for each of the materials as those- given above. 

For properties like those of the boron/epoxy composites of Figure 28, 
the conversion of a small percentage of the axial filaments into transverse 
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(filaments which effectively stiffen the binder should be advantageous, raising ' 

the compressive strength from 2.75 (392 ksi) to 3.8 -S.— (550 ksi), 

m 2 m 2 

approximately (cee Figure 28). If the properties are more like those of the 
.glass/epoxy composites of Figure 29> however, the same gain is not to be 
expected (see Figure 29). If glass is added to boron/epoxy, on the other hand, 
as in Figure 30, substantial gains are to be anticipated. 

Because of the gains found in the calculations for Figure 28 > a series of 
uni -directional boron/ epoxy compression specimens having various percentages 
of chopped-glass filaments added to the binder were fabricated and tested. The 
results are shown as the points in Figure 31 , with the large and small circles 
representing nominal glas s -filament lengths of 0. 16 cm. (l/l6.in. ) and 0. 08 cm. 
(1/32 in.), respectively. Despite the large scatter, strengthening appeared to 
be achieved in two cases. 


Both the scatter in the test data, and the sensitivity of strength of the 
boron/ epoxy composites of Figure 28 to the small amounts of binder reinforce- 
ment point to the critical nature of the role played by the binder in the composites 
for the development of high compressive stresses-. The marked difference be- 
tween the curves of Figures 28 and 29 suggests that this criticality is substan- 
tially greater for boron than for glass reinforcement. 

To explore the sensitivity of boron composites to binder properties 
somewhat further, the calculations represented by the curves of Figure 31were 
also made. These curves show the effects of changes in binder strengths 



and binder-filament-reinforcement strengths on the compressive 
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strength of the three -constituent composite when the "binder strength" mode is 
critical. The compressive strength is found to depend most critically on the 
binder strength at the low ratios of binder reinforcement. 
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STUDIES OF MATERIA LS-PROPERTIES TESTS 

The testing of filamentary composites to determine their mechanical 

•i ' 

properties has proved to be more difficult than, the testing of homogeneous 
materials like metals. The. difficulties take several forms; they may derive 
from discontinuities encountered in specimens designed to provide a diminished 
cross-section, as is often done to avoid failures at points of load introduc- 
tion; they may be associated with the shear weakness of the binder leading 
to undesired modes of failure; or they may arise because a type of test 
known to be inadequate is used for economy or other reasons. In this section 
some of the aspects Of mechanical-property testing of composites are 
examined both analytically and experimentally to help establish'guidelines for 
improved techniques. 

The NOL Ring Split-Dee Tensile Test 

The N. O. L. ring "split-dee" tensile test, shown schematically in 
Figure 32, has the advantages that the specimen may be readily fabricated by 
winding, and the test may be performed in a conventional universal testing 
machine without special-fixtures (other than the split dees themselves). The 
test has the disadvantages, however, that (1) no test section is available at 
which strain gages maybe mounted to measure the stress-strain properties, 
and (2) more seriously, the test introduces substantialbending moments in 
the ring where the split occurs between the two dees. These bending moments 
extend above- and below the, split, as shown in Figure 33. Because the mem- 
brane stress in the ring causes it to increase in diameter as the load increases. 
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the jfing tends to pull away from the corners of the dee, and the maximum 
bending moment is that associated with the load times the deflection to the 
center of the ring cross-section near the corners of the dees. The magnitude 
of the bending may be found from the following equations (developed in Appendix B). 



may be evaluated by the trial and error solution of the equations 


where 




8 and 8 


G 


12 


bending moment at center line of ring, at split in dees 
Young's modulus in direction of filaments 
moment of inertia of ring cross-section 
initial inside radius of ring 
ring thickness 

load acting in each half of ring (one-half of applied load) 
ring cross-sectional area 
unsupported semi-span between dees 

slope at point of departure of ring from dee, measured 
relative to direction of load application. 

■ shear modulus along filaments 


Evaluation of equation (1) for a typical E-glass reinforced epoxy ring 


for R =7.30 cm. (2. ,875 in.) yields the curve of maximum bending stress 
o 

vs. ring thickness given in Figure 34. Obviously these bending stresses are 
of sufficient magnitude to raise questions about the engineering merit of the 
tensile strength values measured as the ^ A stress at the maximum load on the 


ring, 

(Curiously, the magnitude of the bending moment is relatively insensi- 

\ 

tive to the composite properties. For example, the stress increment due to 

bending in an hypothetical boron/epoxy ring with E^ = 414 GN/, m 2 ('60, 000, 000 

_ GN 

psi) is essentially the same as that in a glass-epoxy ring with E^ = 51.75 — ^ 

m 


95 



(7, 500, 000 psi) at the same load. The increased bending stiffness of the boron/ 

.epoxy is compensated for by the decreased diametral expansion of the ring 

v ’ 

under load, so that the bending stresses remain^ssentially constant as 
varies. ) 

Variations of the Split- Dee Test 

The high bending stresses encountered in the split-dee test are rather 
disappointing in view of the many attractive features of the test. Accordingly 
the question naturally arises as to whether minor changes may be made in 
specimen or fixture design which will reduce or eliminate the bending. Some 
possible variations of this nature are considered below. 

The "racetrack" specimen - The provision of a straightaway section adjacent 
to the split in the dees is a first logical step toward the improvement of the 
N. O. L. ring split-dee tensile test. Such a specimen is shown schematically 
in Figure 35. 

The analysis of the bending of the critical section of the racetrack is 
similar to that of the ring, and is described in the following equations:: 
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again with 



'where s = length of straightaway, and the other symbols are as before. 

Results from the analysis of this racetrack specimen compared to those 

/ \ 

from the N. 6. L. ring are compared in Figure 36. Figure 36 shows'that while 
bending is not eliminated'by the straightaway, it is substantially reduced, - 
even by a relatively short straight section. For example, just" a 1 cm (0.4 in.) 
straightaway reduces the bending moment, to less than one-half that of the 
circular ring. 
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The mechanics of the redaction of maximum bending moment are 
suggested in the (exaggerated) sketch of Figure 37. Under load the mid-point 
of the straightaway deflects inward toward the center of the track, so that 
even though the curved portions pull away from the dees the eccentricity of 
■ the tensile load at mid- span is small and hence the moment associated with 
it is small, - i; e. the maximum moment is still that near the corners of 
the dees, not that at mid- span. The fact that the moment in the straight 
portion is small, however, is advantageous for the use of strain gages in.that 
region. 

Experimental confirmation of a qualitative nature of the foregoing 
analyses of split-dee tests were obtained by photoelastic tests (Figure 33 
and 38), and by strain measurements on an enlarged, aluminurm-alloy ring 
(Figure 39). Even with the enlarged ring the stress gradients were so steep 
around the circumlerence that the gages could not be located accurately enough 
to make a quantitative check of the analyses. Indeed, the strain gaging 
problems pointed up the merits of some such approach as the following for 
further improving the split-dee test. 
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Compression Tests of N.O.L. Rings 


A simple compression test for an N.O.L. ring -type specimen (analogous 
to the split-dee test for tension) would be useful for the evaluation of the com- 
pressive properties of filament -wound composites. Efforts to develop such a 
test under this contract were not successful. The specimen tried is shown 
schematically in Figure 40, consisting of two short segments cut from an N.O.L. 
ring and mounted back to back as shown. 

• Application of tangential end loads to the segments of ring in Figure 40 and 
pressure -like forces from the supporting, hour-glass -like core ideally would 
reproduce the compression induced in an entire ring by e'xternal pressure. 

That is, moment equilibrium at any station along the segment is represented 
by the equation . _ 



where 


P * = 1/2 total load 

R. a ring radius 

'a 1/2 angle included by segment 
' a intensity of normal load on ring at station 
that is 

As might be expected, this equation is compatible with. a simple pressure 

• 

load. Thus, if ■*? 
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Tests of specimens fabricated as sketched in Figure 40,. however, were not 
successful in maintaining the compatibility so readily expressible in equation 
form. Attempts were made to replace the resin filler of Figure 40 with aluminum- 
leading to a nice test of the aluminum (Fig. 41); attempts were made to test 
specimens of larger radius of curvature in this fashion, and a maximum failure ■ 

stress of 59 (gg ksi) was so achieved in a glass/epoxy composite. Despite 

mr 

the-fact that failure in this case appeared to be by shearing at the quarter point 

\ 

of the specimen with no evidence of overall column instability (Fig. 42), the 
stress is less than can be achieved with well collimated, straight compression 
specimens and is probably not representative of the material strength. 

Some evidence of "brooming" of the ends of the specimens was apparent 
even in this best test. A number of further tests of like specimens with ends 
cast in cerrobend were made to try to avoid such end failures v Stresses in ,no 
case were higher than the foregoing. value . 1 v ' 

Initially straight specimens appear to be required for compression testing. 
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.Low-Melting Al loy Casting Fixture of Compression Tests 

The use of a fixture in which the specimen ends can be cast in a low- 
melting point alloy has proved successful in preventing premature "brooming" 
type failure in tests of 'straight compression bars. Such a fixture is shown in 
Figure 42. The fixture incorporates end plates having holes, - of the same 
shape as the specimen cross-section, - but approximately 0.6 cm. (1/4 in.) 
oversize so that . 3 cm. of the low-melting alloy can be cast all around between 
the hole wall and the specimen. 

Conventional 3 to 1 aspect ratio, uni -directionally reinforced compres- 
sion specimens have been tested in this type of fixture with consistently high 
stresses at failure. Consistency and modes of failure both suggest that these 
high stresses are representative of the compressive properties of the material. 
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.Elastic Moduli 


The desire to measure the elastic constants of a uniaxial fibrous com- 
posite by using specimens which represent the filament .winding fabrication 
process has led to the use of thin-walled circumferentially wound tubes 
(e. g. Ref 29). Axial load and internal or external pressure tests of such tubes 
can be used to measure four of the five independent effective elastic con- 
stants of the composite material. (See Ref. II for a discussion of independent 
effective moduli of composites. ) The use of these thin cylindrical tubes to 
measure the elastic moduli of a composite laminate appears to be a most 
suitable approach. However, it is not without 'problems; the first of which 
is that of definition of the moduli. When a laminate is to be used as a plate 
or shell structure, the desired elastic relations may be written with respect 
to the principal geometric axes, denoted as the 1 and 2 directions, as; 



where: e are the three middle surface strain components 

X are the three curvature components 
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N are the three stress resultants 
M are the three moment resultants. • 

In this form, for an arbitrary laminate, the C.. matrix is a 6 x 6 array of e- 
lastic constants which must be determined experimentally. Because the lami- 
nate is non-homogeneous there may be coupling between extension and bend- 
ing. The existance of certain non-zero terms in the C matrix complicates 
the relationships among these elastic constants and the so called engineer- 
ing constants: Young's moduli and Poisson's ratios. This is best illustrated 
by considering alternate forms of Equations (8) - (13). These Equations may 
be written as (e. g. Ref. 2): 


N. 

= a 

e. 

+ b.. 

K. 

i 

ij 

J 

.ij 

J 

M. 

= b 

e. 

+ d 

7i. 

l 

ij 

J 

ij 

J 


These Equations are represented by: 



Equation (15)' can be manipulated to yield: 


M 


■ and 


(. 14 ) 


( 15 ) 


( 16 ) 
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(17) 



As an example we now consider the Young's modulus in the 1 direction 
which can be defined as the ratio of average stress to strain and is therefore 
given by: 


N 


t e. 


(18) 


where t is the laminate thickness. 


This modulus can be evaluated from Equation (17) for N = N = M =0, 

2 3 i 


E, = 


1 tA,. 


(19) 


Or it can be evaluated from Equation (16) for N = N = a. = 0 as: 

^ O 1 


E 1 = TT (20) 

1 te u 

The definition of Equation (19) is consistent with the usual practice of having 
all but one of the stress components vanish. However, the restrictions assoc- 
iated with the definition of Equation (20) are representative of the loading con- 
ditions on the shell specimen. That is, a =0, because of symmetry. These 

i ' 

two definitions are not the same unless b . = 0; in which event, it follows that 

U 
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f = g.-. = B = C =0 and A = e 
ij iJ ij ij , ij k J- 

• It is suggested that the rational approach to this problem is to de- 
emphasize the calculation of Young’s moduli,, etc. and to utilize a six by six 
matrix definition of the elastic constants of a laminate. When a particular 
extensional or bending stiffness is required, the definition thereof should be 
explicitly stated. 

A second aspect of the modulus problem is the influence of the form of 
the material anisotropy and the specimen -configuration upon measured values. 
The filament wound cylindrical tube under consideration is a material possess- 
ing cylindrical anisotropy. In the case of a circumferential or other winding 
pattern symmetric with respect to the principal geometric axes, the material 
may be considered to be cylindrically orthotropic. An analysis of the stress 
distribution in such a medium for various applied loads is available in Refer- 
ence 30, which shows that the stress distribution even for simple applied loads 

is non uniform. For example the axial -stress , o , in a cylindrical tube sub- 

z 

lectedto an axial load is not constant unless V = V . This condition is gen- 
J rz 0 Z . 

erally not satisfied in a filament wound shell, and for a circumferentially 

•wound shell, the two Poisson ratios are definitely unequal. In this case, the 

modulus defined by Equation 18 is an average value and can be related to the 

true value by using the expression for the axial stress, a . This result is 

z 

presented on page 254 of Reference 30« However, there is evidently a typo- 
graphical error in Equation (40.4) which is not diemsnionally correct as pre- 
sented. The desired result was therefore rederived for the particular- type 
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of material symmetry -represented by the circumferentially wound tube. The 


extensional stresses and, strains for 

•4 

given by: 



The shear strains and stresses are 


N 

CD 

>- 

= T e z /C 'e 

Y 

= r /G 

zr 

zr zr 

CD 

= T r0 /G r9 


cylindrically orthotropic tube are 


V 



imply related by; 


( 22 ) 


These relations can be further specialized for the circumferentially 
wound tube inasmuch as the material plane normal to the local fiber direction 
can be assumed to be a plane of elastic isotropy. Thus 


E = E 
z r 


V ~ V 

re ze 


G e z = G re 


( 23 ) 
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Equations (21) and (22) are also applicable to filament wound tubes which 
are symmetric with respect to the longitudinal and circumferential direc- 
tions and which have a sufficient number of laminae such that coupling between 
extension and shear is negligible- To take advantage of this wider applicabili- 
ty the analysis will riot utilize^the simplifications offered by Equation (23). 
Since the applied loads are symmetric, there will be no shear stresses or 
strains and the constitutive relations for this problem are those of Equation 
(21). By symmetry, all variables will be functions only of the radial coordi- 
nate. Thus the equilibrium equations simplify to: 


da 


a -a 

^ + -U_§_ =0 


(24) 


which can be satisfied by selecting a stress function, cp, such that 


a 


1 dcp 
r dr 

d 2 cp ' 

a e 2 

dr 


(25) 


The only compatibility equation which is not satisfied automatically 
de 


2 de e 


= o 


(26) 


dr dr \ dr 

.Substitution, of Equations (21) and (25) into (26) yields a fourth order differen- 
tial equation on cp given by: 

4 IV ,3 III 2 II .2 1 „ 

r cp +-2r cp - r gcp + r|3 ip =0 (27) 
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(28) 


, 2 S r 

where g = . — 


»• . s 0 (\-v n y_ 


1-V v 
rz zr 


^0z z6 


E 0 


For p 1 the solution is: 


<p = A r 1_P + Br 1+P + Cr Z + D 


The stresses resulting from this are: 


(29) 


CT = A(l-g)r" (P+1) + B(l+g)r P_1 + 2C 
r 

a. = Ag(g-l) r ' (P+1> + Bg(g+l)r g _1 + 2C 
0 


(30) 


The cylinder is considered.to be subjected to an axial shortening re- 
sulting in'a uniform strain, e = e, and internal and external pressures, p. 

Z ■ v. 1 

and p » respectively. This, train displacement -relations simplify to 
o 


i ( 31 > 


Substitution of Equation (30) into Equation (21) defines the strains, e "'and e Q . 

r 0 

Equation (31) offe'r‘ two 'methBds^of^finding'u, - ' ‘ For ''a 'unique - det'erm'inatibn:' 



Use of the boundary .-conditions: 


108 



a (r=r.) = - p. 
r i 1 


CT J r=r J 

r o 


" P. 


defines the remaining two constants in Equation (30). 


a 

r 




1-C 


2g 


- P 


o 



o 


This yields: 





1-C 

1-C 


P+1 

W 


-(g ii) 

(t) 


1-C 

1-C 


3-1 

2T“ 


r. 

where C = — — 

, T 

O 



P 

z 


[ g -i 

r ml 

jiW-f:) 

1 + C 

1-C 2P 



ri-c 3 - 1- 

i 

Li.c^J 

i 


a 


z 


= E e + v 
z 


0 z 


a + y 

o rz 


cr 

r 


( 33 ) 


(34) 


(35) 


(36) 
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■Axial Compression Test 

It is now possible t,o compare the results obtained by the tube tests with 
the actual elastic moduli. First, consider .the tube in axial compression. The 
measured modulus is given by: 

E = ' (37) 

z e 

a is the average axial stress (total load divided by 
z 

tube cross-sectional area) 


Thus: 


_L 

A 


a dA 


O' 

z 




rdr 


(38) 


Substituting Equations (36), (34) and (35) with p. = P q = 0, into Equation (3? , 
and the- result into Equation (37) yields, after much manipulation: 


E z- E z eV + i) 

1 213 

(1+C 2g )(l+C 2 )-4C P+1 

i<39) 

(1 - y rz y zr> <8 2 -l) 2 

;i + g 2 

_ (i-c 28 

(1-C 2 ) ■ 


v 


For a thin shell we find by taking the limit^of Equation (39) as C-*l that 

E = E 
z z 

Thus the thin-walled circumferentially wound tube appears to be a suitable test 
technique for measuring the transverse Young's modulus. ’ The magnitude of the 
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error associated with this technique is indicated by computing the quantity Y 
in the following equation as a function of p : 


E -E 
z z 

E 


z 


l-v v 

rz zr 


(40) 


where the correction factor, ¥, is determined by comparison of Equations (39) 
t ' 

and (40). 

• The factor ¥ was evaluated for a radius ratio of C = 0. 9 which is a rela- 
tively'thick shell for this purpose. Correction factors, ¥, for smaller radius 
ratios, C, will be smaller. It was found that’: 


¥ < 0. 01 2,<PS6 

¥<0.1 2 < (3 < 10 

It seems safe to conclude that the circumferentially wound tube in axial 
compression is a valid test for the Young's modulus transverse to the fibers of 
a uniaxial composite. 
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Internal Pressure Test 


An internal pressure test has been used (e. g. Ref. 29) to determine the 

' V 1 

Young's modulus in the fiber direction and the Poisson's ratio. Strains are 

j 

Pleasured on the inner and outer surface. For this loading we may examine 

;the equations for strain by substituting Equations (34) - (36) into Equation (21). 

For zero longitudinal strain and zero external pressure this yields: 
i . 


-(PH) 





201 


9r 

'ZG 



(41) 


The measured surface .strains will then be' given by: 


® 0 (r i> = 


P,- 


1-c 


20 


1+c 


20j 


9 r 


(l-c 2P ) 


(42) 


0+1 


^e^o 5 = 20 le P i 29 

1-C 


(43) 


As in the preceding section, we choose to examine the error associated 
with the use of equations developed for isotropic materials. The alternative 
is to write the simultaneous equations relating all the measured strains to all 
•the elastic moduli and solving for the latter. However, this alternative pro- 
cedure is complex and perhaps not possible as tliere is no apparent direct 

measurement to define the transverse shear stiffness, G . Thus, we con- 

, zr 

sider 
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(44) 



where 


For equilibrium: 


Eq is the effective circumferential thickness for 
' zero axial strain. 

is average circumferential stress. 


3 


0 



(45) 


The stiffness defined by Equation (44) can be evaluated by the use of e. g. (45) 
and either Equation (42) or (43). We denote the results Eg and 1° respec- 
tively.. Treating the thin shell we write: 


c = 1-6 • (46) 

where 6 is the shell thickness to radius ratio. The results are: 


E , 


r 


[-41 


(47) 


E 1 _ _L 

•E* ft c 


1 -- 


L+_l5l 

2 G S 


rs 8 j j 


(48) 


Equation (47) shows that the use of the external surface strain measurement 
yields a result for the extensional stiffness which has an error measured by 
half the thickness to radius ratio, e. g. 5% error for a radius to thickness 
ratio of 10. On the other hand, the inner surface measurement involves many 
other material moduli and is therefore unsuitable for simple use in the circum- 
ferential modulus determination. 

Note also'that there is a significant difference between internal and ex- 
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ternal strains. Thus: 


e e^ r i* " e s^V 




ier$r»* 


G r Z s e/ 1 “ (p+1)6 


Reference to Section 43 of Reference 30 'shows that the stress distribu- 
tion in a cylindrically anisotropic hollow shell for which there exists at each 
point a plane of elastic symmetry normal to the axis is obtained in the same 
way as in an isotropic rod. That is, ^ is the only non-vanishing stress 
component and 

M 

a 6z~ I - r 
o 

The circumferentially wound circular shell satisfies these symmetry require- 


ments and therefore provides an exact measure of the in-plane sheaf stiffness. 
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CONCLUDING REMARKS 

• -The application *of existing analyses of strength and stiffness of fibrous 
composites to the analysis and design of composite structures has been treat- 
ed . It has been shown that, although the present understanding of failure 
mechanisms is incomplete, it is possible to construct a rational set of fail- 
ure criteria for fibrous composite laminates. Additional studies of the fail- 
ure mechanisms of a uniaxial composite have been performed by using limit 
analysis techniques. ‘Also the effect of stress concentrations upon tensile 
strength has been explored and initial 'studies were made of the complex 
moduli of fibrous composites having viscoelastic matrices. 

The design criteria were used in the evaluation of columns, plates and 
panels in which uniaxial composites for high strength were combined with 
isotropic metals or composites -for high stiffness. Designs offering substantial 
improvement over those made of one material only were demonstrated. Im- 
proved composite materials were achieved through the use of a third phase as 
a local reinforcement of the matrix. Also approximately isotropic three- 
dimensional materials were designed. 

Test methods were analyzed leading to a verification of the suitability of 
thin-walled tubes for modulus measurements and the unsuitability of the NOL 
ring for tensile strength measurements. Improvements in methods for measure 
ment of tensile and compressive strength were demonstrated. 

The studies emphasize the feasibility of using present methods for the 
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mechanical analysis of composites in preliminary design studies. They also 
|ndicate that our understanding of composite failure mechanisms is a continual- 
ly evolving one requiring additional theoretical and experimental study. 
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APPENDIX A - Limit Analysis For Composites With Elastic -Brittle Fibers 

; Five types of surface loadings are considered: 

1. • Tractions equivalent to' a uniform shear stress -T^ are applied on 
the entire boundary surface. 

Use 'will be made of the theorems of limit analysis which will not be 
stated here. For lower bound construction, a. uniform stress field 

0 r 0 

o 

T - 0 0 

o 

0 0 0 _ 

is chosen as a statically admissible stress field where T is such that Von 

o 

Mises ’yi e l.cl criterion of the matrix is not violated. Then it can easily be 
shown that 



r. . 

..ij 


where k^ is the yield stress of the matrix in simple shear. 

For upper bound construction, a kinematically admissible velocity field 
is chosen as follows: 

(a) In the region of the composite specimen not occupied by composite 
cylinders and on the boundary of the composite cylinders. 


U 1 = °> 


u = 3 

2 '0 


u 3 = 0 


( 2 ) 
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where Y is any real number’. ’"' 

0 

(b) In any composite cylinder, the velocity field vi is the elastic dis- 
placement solution to the displacement boundary value problem with displace- 
ment boundary conditions (Eq. 2) prescribed as formulated in Appendix 2 of 
(Ref. 11) with the modification that the fibrous core is rigid. 

Ah application of the upper bound theorem gives an upper bound (T ^) 

• 1C U 

as a function of P and v^ the latter being the volume fraction of all the composite 
cylinders embedded in the composite specimen: 


r u) u 3 V Vi 'V 11 


(3) 


where 


-4-f f 

1 J 5 J 0 R R 


cos 0 d0dR 


In the case of '’random array' 1 (Ref. 11) for which v = 0 and the fiber 

r ^ 

2 - 

volume fraction of the composite specimen v^ = P , Equation (3) becomes 


12 , 


U 


s = I > 


(4) 


2. Tractions equivalent to a uniform shear stress r ^ ar S applied on 
the surface of the specimen. 

The condition of plane strain is assumed. Then, the Von Mises 1 yield 

criterion for the matrix reduces to a simple form 
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( 5 ) 


( T 22 ' T 3 3 ) 2 + , *h\ 2 4U b 

r 

For lower bound construction, a uniform stress field 


0 0 0 



pro 

o 


is chosen as a statically admissible stress field where r is such that Eq. (5) 
is nowhere violated. Then it follows that the lower bound for the limit load 



which is independent of.fiber volume fraction. 

.For upper bound construction, the same principle is used here. A 
kinematically admissible velocity field is chosen as follows: 

(a) In the region of the composite specimen not occupied by composite 
cylinders and on the boundary of the composite cylinders. 



(b) In any composite cylinder, the velocity field u is the elastic dis- 
placement solution to the displacement boundary value problem with displace- 
ment boundary conditions (Eq.' (6) ) prescribed as formulated in Appendix 1 
of Ref. 11 with an additional condition that the fibrous core is rigid and the 


1 22 ' 



binder shell is incompressible. 

An application of the upper bound theorem gives ah upper bound 
as a function of •$ and v : 




= k b +k bV I 2- 1) 


(7) 


where 


tt(1- 


3 2 ) 3 


■1 TT 


K) 


i)i 2 (p, R) cos 0 d0d R 


and 


P 4 

R 4 u 


r y 21 

2 

F . 4 2 ^ H 

L R. 

+ 

-( 4P 2 +P+ 1» + 3f?R 2 


<L(P, H) = % 
2 ' R 2 


(pV + l) -|-(P 2 +1 )^ 


R 2 -JL 


(4g 4 +p 2 +l) 

2 


2 2 

+ 3p R 


3. Tractions equivalent to uniform axial tension t on composite • 
specimen under the conditions of plane strain. 

Using the same principle, the lower bound for the limit load is 



= 2k b 


For upper bound construction, a kinematically adraiss ible' velocity field 
is chosen which is obtained from the one constructed for the preceeding case 
through an orthogonal transformation such that 
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u = d u 1 with d = 0 

1 ij j ij 


0 ■ 0 

1 1 

a y 2 


0 



_i_ 

ifZ 


where ul are the velocity components used in the preceeding case referred to 
an x' -system. Then, after some manipulation, it turns out that 

l\ 


22 


= 2 t. 


23/ 0 


Another kinematically admissible velocity field was constructed: 
a. In the region of the composite specimen not occupied by composite 
cylinders, and on the boundary of the composite cylinders, 


V 0 - U 2 = TV U 3=-T X 3 


or in polar coordinates. 


' Y q Y 

u = 0, u = — r , cos 20, u. = - — r sin20 
1 r 2 b 8 2 b 


b. In any composite cylinder, the velocity field in polar coordinates 
is assumed to have the following form: 

Y 


u = — p(r) cos 2 0 
. r 2 


for r £ r £ r ■ 
f b 


2 - 


Q (r) sin2 8 
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where p (r) and ft(r) are arbitrary functions and 


p(r b ) = n( r b ) = r b 


( 8 ) 


In addition, we assume that p(r^) = fl(r^) = 0: Becuase of the assumed 

incompressibility of velocity field, p(r) and Cl(r) satisfy the following equation: 

' d P(r) 


d r 


2 0 (r) + p(r) 


0. 


( 9 ) 


2 3 

Furthermore, let p (r) = A^r-r^) + 

a '.ilils 

■ 1 ' (r - . ) 2 

b f' 


where the constants 


and 


A 2 = ' 


r f + r b ^ 
(? b' r f )3 


have been determined from (8) and (9). 

Now the velocity field of the entire specimen is constructed and the 
upper bound theorem is applied to obtain a new in terms of fiber- 

volume fraction for the case of "random array". Incidentally, this same vel- 
ocity field, after an orthogonal transformation, can also b‘e used to obtain a 


Wt 


new{r.l . Unfortunately, numerical calculations show that for all fiber- 

7q n 


volume fractions between 0 and 1, bothfrll and (r I obtained here' are 

\ 22 /u ■ V23^ ■ • 

slightly higher than their corresponding values obtained previously. 
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4. Tractions equivalent to biaxial uniform tension t and r 
i' 22 '33 

(j_ 7 s T 33 ) on the composite specimen. 

For definiteness assume r 22 > r 33 - Then following the ideas used in ' 
Cjase 3, it can be shown that 

( ■ ( T 22^ ' ( T 33^ = 2k b 


5. The type of surface loading considered here is for r. . in the follow - 


Lng. form; 


T 12 T 22 


where the constant stress components 7^ and r are related in the following 


T 12 = “ T 22 


with Cc ^ 0. This amounts to a p r oport ional loading of combined uniform shear 

stresses Tj 2 and uniform tensile, stresses t ’ on the boundary surface S of 

\ 

he specimen. 

Since both j ^ an< ^ ^2 are assume< ^ finit'e, it is obvious that a= 0 
:orresponds to the case where only uniform tensile stresses a re present. 
)n the other hand, a-> ■» corresponds to the case where the specimen is sub- 


126 



iected to uniform shear stresses r 

12 

In lower bound construction, a uniform stress field of the form (2) 
supplemented by (3) is chosen as a statically admissible stress field to obtain 
lower bounds for the limit loads and , respectively. It turns out that! 


// L»\ 

' V 22 ) L 


for as= 0. 


From Equation (12) 'it 'is obvious that for the special case where a = 0 
(i. e. for uniform transverse stress r • acting alone), 


In the case where a-» ", corresponding to the case where the specimen 


is subjected to the uniform in-plane shear stresses t a l°ne, 


lim \12 
a-» " k 


/ L\ 

& = 1 
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; For upper bound construction, a kinematically admissible velocity 

field u is chosen to be' a linear combination of the two kinematically ad- 

•■9 u 

u ■" u 

(a) In the region of the composite specimen not occupied by composite 
cylinders and' on the boundary of the composite cylinders: 


missible velocity fields used to obtain and (r,V) respectively: 

< ' 1 2 4t ' TT 


U) X (2) 
u .+ u 


(14) 


where 


and 


u (1) = ( 0 , y x . o) 
-> 11 


( 2 ) 


Y v 

/ 2 2 \ 

V 0 ’ T x 2’ - T x 3 ) 


with . 


Y 


1- 


U) y 

2 


(15) 


where UJ is any real number. 

/ >. « 

(b) In any composite cylinder, the velocity field is the elastic displace- 
ment solution to the displacement boundary conditions ( 14 together with (15) 
prescribed with an additional condition that the fibrous core is rigid and the 
binder material is incompressible. In fact, the solution to this displacement 
boundary value problem can be obtained from the associated solutions to the 


individual problems connected with obtaining 



by the prin-. 


ciple of superposition. ' 

Then an application of the upper bound theorem gives,, for' the case of 
constant g throughout the specimen 
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( 16 ) 


T zzj u 


•where 


,1 „2tt 


3 tt(X-P 2 > 3 


n (U), p, R) = 4, 4 

1 • . ( R 


R^JQ^ui, p, R) + £!.,(«>, p, R) cos 0 + (p, R) cos 20d0dR 

(p 4 +p 2 +l) (p 2 +l) 


„2l 2 


2 2 
3p- R 


2 2 
+ u) (-1-P ) 


(4g 4 4 g 2 + i) 
2 


Q. (IB, P, R) = 2(l) 2 (l-p 2 ) 4 -^ 
2 R 


and 


• Cl (P, R) = 8^ 

3 R 2 


(p 4 +p 2 +l) - — (P 2 +l)- 


2-i r 


2 2 (4P 4 + P 2 ' + 1) 


3 P R - 


' ' 2 J 


In the case o£ "random array", (16) reduces to 



k 1 

b j + a U) 
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Since U) in (15) is arbitrary, ■ the lowest upper bound among the class 


o| upper bounds in (17) will be obtained by minimizing the right-hand side of 
(.1-7) with respect to u>.’ 

| Thus, 



( 18 ) 


and 



kal 
b 3 

“ +au) 
2 


will be chosen as the upper bounds for 


22 


and r 


L 

12 ’ 


respectively. 


Numerical calculation is performed to obtain (r^j and from 

(18) for. different values of g and a. It is interesting to note that in the num- 

I 

erical calculation, for any given g, U) which minimizes the right hand side of 
Equation (18) is a monotonir; increasing function of a but U) i a (except when 
a = 0, then' u) '.= a = 0) . 
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APPENDIX B - Derivation of Equations for N. Q. h. Ring Split- Dee Test Analysis 
The assumptions are made that the N. O. D. ring exactly fits around. 

. V 

. the two dees at zero load and that there is no friction between dees and ring, so 
t|iat the portion of the ring that remains in contact with the dees is subjected 

i 

to essentially a radial pressure which induces a hoop tensile force P in the 
ring. This hoop tension produces a circumferential expansion of the ring, 
equivalent to: 

8 

where 

Young's modulus in hoop direction of ring - 

A cross-sectional area of ring 

R initial inside radius of ring 

o 

t ring thickness 

This circumferential expansion is converted to a straight-line separation 
Ap between the dees, so that 



very nearly, and the problem is that of the analysis of abeam at the split in the 
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dees as sketched below: 



That is,' the beam is tangent to the dee at some point at a distance 

from the center line of the split between the Dees. At the centerline the beam 

is subjected to the tensile force P and a moment . At the point of tangency 

the force P. is also acting and a fictitious moment is hypothesized of the 

magnitude required to form the ring to .its initial radius of cur-vature R + 4r . 

. _ o Z 

So 



In other words the analysis considers the beam as initially straight when un- 
stressed. The desired bending moment at the center line M will thus be found 

I XJ 

by subtraction of the fictitious moment M from M, . ’ * 

o 

Selecting x-y coordinates as shown, we note\that the moment M at any 


section of the' beam is . 



or 



( 4 ) 


Let 



By substitution (4) becomes 


.(5) 



The solution of (6) is 




■ K ) 


fax 


(7> 


Accordingly the slope @ at the point of tangency of the beam to the dees is 


^ sink 




7 


The deflection £ is similarly found as 
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Combining (8) and (9) 



Because the beam is tangent to the dee at x = A .the angle (Q subtended by 
the arc of the dee between the split and the point of contact may be found by 


. trigonometry -as 


/ W 


7 l 


Finally, from (3) for y = 0 


Or by substitution of (9) 
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Shear deflections may be included in the analysis by re-writing equation (7) as 



For a rectangular section.n = 1.2, and (15) becomes 



Analysis to include the effect of the straightaway on the racetrack specimen 
follows the same form as above, leading directly to the equations given in the 
text. 
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/- Direction Ft /aments 



Generalized equations for compliances of composites having uni- 
directional reinforcing filaifients in the 1-direction (from Ref. 1). 

136 



Table 2. 


h Direction Filaments' 

In ierrelathmhips among constants for 
transverse isotropy • 








A 


7 



S 



Relationships among elastic constants and compliances for uni-direc- 
tional reinforcement for use in evaluations of the various « 's in 
Table 1 (from Ref. 1). 
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Table 3. - Equations for the transverse effectiveness factors ($) in terms of 

the elastic constants for uni-directional reinforcement (equations of 
Table 2 solved in terms of gfor 0 values equal to unity). 
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£ i ) fi O'^i) / \ ■ 

J &/'-%) „ . ^( h x i „ L ) 

I Eh(>-‘i.) L ) 

f j o**f 1 Xi-*%) t3 ’ T ’ ’ O'nX'r^dP-’ipr . 




to 



Table 4. - Generalized equations for the compliances of three-dimensionally 
reinforced composites having three orthogonal planes .of symmetry. 
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•Oa E>{_ ^b ^b / ) 

+ fi+ % x 
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Table 4. - (Cont. ) 
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A 





A 9 " G i pX- + 1 G fX, ' 




cos*Jl/co$ty-j//> A * co$*<{> costy 


sin^d 

!*%] 


hJ 
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Table 4. - (Cont. ) 
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STRESS, KSI 



Figure 1. 


Calculated Stress -strain Curves for E -Glass and Epoxy Composite 
Laminates. 


N 401-923 


Figure Z 


Yield Strength of a Symmetric Bi -axial Composite Laminate for 
Failure Mode's Involving Each of the Principal Lamina Stresses. 
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Figure 4. Bounds on the Limit Load for In-plane Shear, 
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Figure 5. Bounds on the Limit Load for Transverse Tension. 
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NUMBER OF BREAKS 



Figure 9. Number of Breaks as a Function of Applied Load (Ref. 4, Series B). 


NUMBER OF BREAKS 



Figure 10. Number of Breaks as a Function of Applied Load (Ref. 4, Series C). 



EXPECTED NO. OF FRACTURE GR 



EXPECTED NO. OF FRACTURE -GROUPS 



Figure 12. Expected Number of Groups of Fractures as a Function of Applied Load (Ref. 4, Series C) 


NUMBER OF MULTIPLE- FRACTURE GROUPS 



FIBER STRESS (KSI) 

Figure 13. Observed and Predicted Multiple -Fracture Groups {Ref. 4, Series B). 




Figure 14. Observed and Predicted Multiple -Fracture Groups (Ref. 4, Series C). 
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Figure 15. Static and Dynamic Values of an ^ ^3/1 aS a function Applied Load (Ref. 4, Series B). 
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Figure 17. Variation of E with Applied Stress for Various Composite Lengths (Ref. 4, Series B). 
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Figure 18. ■ .Typical Variations of Transverse Effectivenesses of Filamentary 
Reinforcement with Volume Fraction, as Calculated for 'an E- 
Glass/Epoxy Composite. 
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Figure 19. Efficiencies of Round, Unreinforced Tube - Columns of 7075 ~T6 Aluminum Alloy and 
Beryllium. 
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Figure 24. Efficiencies of Beryllium Scalloped Triangle Thin-Walled Tube Columns Reinforced at 

Apexes wtth 0 Boron/Epoxy Stiffeners, and Comparison with Minimum- Weight R einf orced 
Beryllium Round Tubes. 
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Figure 26. 


Comparative Efficiencies of Wing Box Beams having ±45° Boron Reinforced Epoxy Compres- 
sion Skins and Beryllium Skins on 0° Reinforced Boron/Epoxy Z -Section Stiffeners. 





f 

A s 

fl'* 

A* s 

V* 

Af 

0.1 

1.39? 

144/- ■ 

1383 

1375 

.4 

£.198 

£234 \ 

£.&§£ 

LUO 

.6 

0.26& 

3.273 

3.034 

3.093 

■ M 

4.077 

<4. $87 

4.525 

4.779 


Figure 27. Calculated Compliances Ap. A^, and A^, Caltropic Reinforcement 

in the Directions of the Orthotropic axis 1-, 2-, 3- Res-pectively, for 
Equal and Unequal Assumed Transverse Effectivenesses /?. (E^ = 10, 

E b = if' V f = °-‘ 2 - ^ = 0-35). 
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Figure 28. Calculated Compressive Strengths for Three Indicated Failure Modes 
for a Boron/Epoxy Composite of Nominal Constituent Properties with 
Constant Total Reinforcement Volume Fraction of 50% but Varying 
Proportions of the Reinforcement Uni -directional and Randomly Dis- 
persed in the Binder. 
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Figure 29. Calculated Compressive Strengths for Three Indicated Failure Modes 
for Glass /Epoxy Composites of Nominal Constituent Properties .with 
Constant Total Reinforcement Volume Fraction of 50%, but Varying 
Proportions of the Reinforcement Uni-directional and Randomly Dis- 
pensed in the Binder. • 
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Figure 30. Calculated Compressive Strengths for Glass -Boron/Epoxy Com- 
posites of Nominal, Constituent Properties with Constant Uni-Direc- 
tional Boron Reinforcement Volume Fraction' V£ = 0. 5 and Varying 
Quantities of Randomly Dispersed Glass Filaments in the Binder. 
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Figure 31. Experimental Results for Tests of Glass -Boron/Epoxy Composites 

having Uni-directional Boron Filamentary Reinforcement of Nominally 
50% by Volume, and Comparison with Calculations for Various 
Glass-Filament and Binder Strengths. 







Figure 32. 


Schematic of N. O. L. Ring, "Split-Dee" Tensile Test. 
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Figure 33. Photoelastic Study of Stresses in 
Ring "Split-Dee" Tensile Test. 


the Vicinity of the Split between the Dees in an N.O.L. 
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Figure 35. "Race-Track", Filament-Wound Tensile Specimen. 
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Ratios of Maximum Bending Moments Induced in Race -Track 
and Circular 0. 15 cm. [0. 06 in.) Thick Split-Dee Specimens of 
Glass/Epoxy at 0. 69 (i 00 ksi) Axial Stress. 






Figure 37. Schematic Representative of Mechanics of Deflection of "Race- 
Track" Specimen. 
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Figure 38. Photoelastic Study of Stresses in the Straightaway of a "Race- 
Track’ 1 Specimen. 
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Figure 41. Aluminum Insert from Specimen Similar to that of Figure 40 after Test. 
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Figure 42. N. O. L». Ring Segment Compression Test Specimen Showing Shear Failure. 
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